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1. INTRODUCTION 


This document presents the final technical report on work performed 
at Virginia Polytechnic Institute and State University under the ^ASA 
contract NAS9-13091. This contract provided for the development of a de- 
tailed and accurate numerical model of the NASA-DELCO Electromechanical 
Actuator or EMA, References (1), (2), (3), (4). 

For convenience, the final report has been divided into two separate 
volumes: Volume I and Volume II Volume I deals with the development of 

a numerical model for the EMA. This topic will be expanded upon in 
Reference (1). The first volume also contains a comparison between data 
obtained from the numerical model and data obtained from tests performed 
on the actual hardware. Volume II describes the function and use of the 
FORTRAN COMPUTER PROGRAM which implements the numerical EMA model. 

The EMA model developed here is based upon the descriptions given in 
References (2), (3), and (4). Background information and descriptions of 
the NASA-DELCO EMA are presented next. 

1.1 THE NASA-DELCO EMA 

At the present time, the positioning of flight control surfaces such 
as flaps, elevens, and rudders, Figure (1.1-1), is accomplished almost ex- 
clusively by hydraulic actuators. The overwhelming preference for the 
hydraulic actuator is due to the poor past performance of the EI-IA when 
compared with its hydraulic counterpart. This poor performance was due 
mainly to two factors: (1) The lack of practical, fast and reliable high 

power transistor switches and (2) The lack of a suitable actuator motor. 

Recent advances in magnet technology. Reference (6) , have made it 
possible to reduce the size of both the rotor and stator assemblies through 
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Che use of high energy density samarium-cobalt permanent magnets to re- 
place Che rotor field winding. The reduction in rotor size and weight 
(mass) has resulted in Improved actuator performance because of the reduced 
rotor inertia. At the same time, transistor switches capable of switching 
more chan 10 KW in one microsecond or less became available, References (2), 

(3) and (4). 

As a result of these advances, a technology program was initiated by 
NASA to explore the possibility of replacing existing hydraulic flight con- 
trol surface actuators with electromechanical ones. The justifications 
for this program are based upon the possibility of reducing actuator weight, 
volume, and maintenance requirements as well as Increasing reliability. 

NASA' s interest in Che EMA concept resulted in a contract to the 
DELCO Electronics Division of General Motors; see References (2), (3), and 

(4) . This contract, which has been completed successfully, provided for 
the development, fabrication, and evaluation of a prototype EMA capable of 
performing in aerospace applications. 

The actual function performed by this prototype EMA is illustrated 
in Figure (1.1-2). In this figure, the actuator receives a flap position 
command from an external source (the pilot for instance). The EMA control 
logic translates this position command into a motor torque of a specified 
magnitude and orientation. This commanded torque is produced by a power 
conditioner which drives an eight pole permanent magnet brushless d.c. 
machine. The resulting flap motion is in the direction of the commanded 
flap angle. Once the flap reaches the commanded position, it oscillates 
briefly about this point until the energy of the rotating masses is com- 
pletely transformed into both heat and electrochemical energy. 
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The physical structure of the Delco EKA is given in Figure (1.1-3). 
Notice that this system utilizes four Independent servo loop channels. 

The four Independent channels help to Increase the overall reliability 
because during normal operation two channels are active while the other 
two are Inactive. High reliability Is of great importance In aerospace 
applications since the failure of such components results In reduced or 
total loss of aircraft control. 

The structure of one of these servo loop channels Is shown In 
Figure (1.1-4). At the heart of each servo loop channel there is a power 
conditioner and machine combination, see Figure (1.1-5), which modulates 
the power of the supply batteries Into the proper actuation energy. In 
order to achieve this goal, the power conditioner and machine network are 
operated In any one of three modes by means of transistor switching com- 
mands generated by the control electronics. These modes are: 

a) Motoring in the Positive and Negative Directions, 

b) Regenerative Braking in the Positive and Negative Directions, 
and 

c) Plugging or Braking in the Positive and Negative Directions. 

The control electronics determine the proper operating mode at ary in- 
stant in time based upon the Inputs shown in Figure (1.1-4). These in- 
puts are: 

a) The Flap Position Command • 

b) The Three Bit Output of the Rotor Position Sensor or RPS 

c) The Machine Velocity 

d) The Flap Position 

e) The Average Machine Current (Torque) 

A typical sequence of mode commands generated during a step response, 
such as the one in Figure (1.1-2), is shown in Figure (1.1-6). At time 
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t ■ t the machine Is commanded to motor In the positive direction in 
o 

response to the step position command. The flap responds 
with motion in the positive direction, thereby reducing the distance be- 
tween the actual and commanded flap positions. When the flap reaches the 
commanded position at some time t ■ t^^, certain logical decisions are 
made. If the rotor velocity at this point is greater than a specified 
threshold velocity, then the motor is thrown into the regenerative braking 
mode. References (l)-(4). During this mode, the machine behaves like a 
constant field current salient pole generator. When the magnitude of the 
machine velocity dips below the threshold value, the EMA enters the plug- 
ging 07 braking mode. The machine In this case receives energy from the 
supply battery which decreases the energy of the rotating masses. When 
the rotor comes to rest, the position of the flap will exceed the com- 
mand posltioii. Consequently, the machine is forced to motor in the nega- 
tive direction. This process is repeated until the commanded flap angle 
is reached ar.d the flap is at rest. 

The velocities of the four machines are summed by a differential 
gearbox arrangement as shown in Figure (1.1-3). The physical details of 
this gearbox are pictured in Figure (1.1-7). The advantage of the dif- 
ferential gearbox is the isolation of the rotor inertias of the inactive 
machines from the inertia of the active machines. The output shaft of 
the differential gearbox drives a flap hinge actuator. The flap position 
is fed back continuously to the channel control electronics by position 
transducers attached to the differential gearbox. Figure (1.1-3). 

The rotor velocity is sensed by means of a tachometer attached at 
the back of the machine. This velocity along with the flap position and 
the average machine current (or Torque) is processed by the control elec- 
tronics into a current command for each machine. 
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The current conmand is equivalent to a machine torque command be- 
cause of the linear relationship between the average machine currents 
and the developed electromagnetic machine torque, due to the equivalent 
constant field excitation, References (1), (2), (3), and (4). 

The EMA described here has been tested and analyzed by Delco 
E3.ectronics, References (2), (3), and (4). The results of these tests 
have demonstrated the potential of such actuators for aerospace applica- 
tions. However, there are still many questions which must be answered be- 
fore EMA's become serious contenders against existing hydraulic actuators. 
For example: 

(a) Can the switching control for the power conditioner tran- 
sistors be modified so as to reduce the severity of the 
switching transients? 

(b) How does the timing of the phase commutation affect the 
overall response and reliability of the EMA? 

(c) What are the effects of system parameter changes on the 
overall system efficiency and performance? 

(d) What are the consequences of component failures? 

These and other questions &re difficult to answer experimentally because 
of the high costs Involved in prototype development. Therefore, the only 
other alternative is to simulate the EMA dynamics by means of a mathematical 
model of sufficient detail to allow the accurate calculation of the EMA 
response under a variety of Inputs and configurations. The development and 
verification of such a model is the topic of the remainder of this first 
volume . 
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1.2 PROBLEM DEFINITION AND STATEMENT OF TASKS 

The EMA which was introduced in the last section must be examined 
thoroughly under all possible operating conditions in order to determine 
the expected ranges of the system state variables. This Information is 
vital in the determination of the stresses experienced by the elements of 
the power conditioner and hence its reliability. Furthermore, the present 
design must be carefully evaluated for the purpose of design optimization. 

Such knowledge is very costly to obtain by experimental means. In 
addition, some of this knowli^dge would require tests which might result 
in damage to the system. For these reasons a mathematical model of the 
EMA becomes a necessity. 

The EMA actuator unfortunately contains many nonlinearities which 
greatly complicate the mathematical model. The differential equations 
which constitute this model are therefore nonlinear. The only practical 
solution of such equations is by means of a numerical integration scheme. 
Based upon these observations it becomes evident that a COMPUTER SIMULATION 
of the EMA on an instantaneous time basis is the only practical modeling 
approach. The chosen modeling approach can be outlined as follows: 

STEP 1 . Determine the set of nonlinear differential equations which 
govern the dynamics of the EMA from the underlying physical principles. 

STEP 2. Simulate the control logic by meant of equivalent boolean 
expressions. 

STEP 3 . Divide the solution period into N equal time increments or 
steps. These increments are chosen sufficiently small so that the non- 
linearities can be linearized during each increment. 


STEP 4 . Integrate the differential equations forward by one time 
increment at a time until the entire solution period has been covered. The 
differential equations are updated, if necessary, after each integration 
in order to account for the nonlinearities. 

These steps are covered in greater detail in Chapters 2, 3, and 4 
of this volume. A model of the power conditioner and machine network is 
derived in Chapters 2 and 3. The equations which describe the control loops, 
the rotating masses and the control logic are developed in Chapter 4. 

The specific requirements of the NASA contract, NAS-15091, are sum- 
marized by the statement of •. asks below: 

TASK 1 : Development of a mathematical representation or model for 

the entire EMA system shown in Figure (1.1-2), under one or two machine 
operations. This model must include: 

(a) A simulation of the logical functions performed by the 
low level electroi.ics. 

(b) A simulation of the linear and nonlinear transfer func- 
tions which form the servo loop controller. 

(c) A simulation of the nonlinear power conditioner network 
under all operating configurations. 

(d) A simulation of the electromagnetic interactions inside 
the permanent magnet machines. 

(e) A simulation of the rotating masses. 

TASK 2 : Implementation of the numerical EMA model, developed under 

Task 1, in terms of a FORTRAN computer program. 

TASK 3: Verification of the numerical model by means of test results 

obtained from the actual hardware. 

TASK 4 : Determination of the commutation transients under various 

operating regimes. 
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TASK 5 ; Development of a USER's guide for the FORTRAN program. 

These tasks have been completed successfully: and together they 
form the basis for this report. Tasks 1 through 4 are covered in the 
first volume of this report while Task 5 constitutes Volume II. 

The upcoming chapter Introduces the actuator motor and models 
which describe its performance. 


2. ELECTRICAL MODEL OF THE PERMANENT MAGNET 
BRUSHLESS D.C. MACHINE 

Chapter 2 deals with the four permanent magnet brushless d.c. 
motors which power the Delco EMA. An electrical network model of 
these motors is developed. This model Includes the effects of magne- 
tizing and demagnetizing mmf's of the stator winding. Later, in 
Chapter 3 this model is Incorporated with the power conditioner model 
to form the power conditioner and machine network model. 

This chapter is divided into four sections. The first section 
describes the mechanical and electrical characteristics of this machine. 
The three different modes of operation: Motoring, Regeneration, and 

Plugging, are described in terms of the electromagnetic Interactions 
which take place Inside the machine in Section 2.2. The electrical 
network model of the machine is developed in Section 2.3 which follows. 
In the last section of this chapter, the calculations of the machine 
back emfs and the developed electromagnetic torque are presented. 
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2.1 MACHINE DESCRIPTION 

The electromechanical actuator, introduced earlier in Section 1.1, 
derives its actuation energy from four identical permanent magnet brush- 
less d.c. motors. This type of machine was chosen because of its relia- 
bility, low rotor Inertia, low weight and volume per KW, negligible rotor 
eddy current losses, and excellent torque characteristics. These factors 
are important if this type of actuator is to compete successfully with 
conventional hydraulic systems. These machines owe many of their excel- 
lent performance characteristics to the recent advances in the field of 
high energy product permanent magnets. Such magnets have made it possible 
to reduce the size of the rotor without risk of demagnetizing the magnetic 
pole pieces, see Reference (2). 

Each one of these four machines has a rating of 17 hp at 9000 rpm 
with an output torque of 13.558 nt-m (120 in-lb). These ratings are 
achieved with a permanent magnet rotor of only 0.1335 m (5.25 in) in 
length and 0.A72 m (1.858 in) in diameter, see Reference (2). More 
detailed mechanical and electrical information about these machines is 
presented in sections 2.1.1 and 2.1.2 which follow immediately. 
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2.1.1 Mechanical Description 

A relatively complete description of the Delco permanent magnet 
brushless d.c. motor can be found In Reference (1.1-1). However, for 
the sake of completeness, a brief description will be given In this sec- 
tion. The major components of the permanent magnet brushless d.c. motor 
are: 

a) The Rotor Assembly, 

b) The Stator Assembly, 

c) The electromechanical Brake, ^ 

d) The Rotor Position Sensor, 

e) The Tachometer, 
and 

f) The Cooling System. 

These components are clearly shown In Figure (2.1-1) which gives a cross 
sectional view along the rotor axis. The outside diameter of the motor 
frame shown in this figure is 0.09525 m (3.75 inches). The total weight 
of the motor Including the brake, tachometer, and rotor position sensor 
is only 9.025 kg (19.9 lbs). 

The rotor assembly consists of 48 blocks of magnetic samarium co- 
balt material which has an energy product of 18 million gauss oersteds. 
These blocks are bonded, six In a row, onto the faces of an octagonal 
rotor shaft. Figure (2.1-2), creating an eight pole magnetic 
structure. Figure (2.1-3). Brass and disks are bonded to this assembly 
to reduce windage losses and to provide removable material for dynamic 
balancing, see Reference (2). The rotor assembly is then banded by 
a fiberglass material. Figure (2.1-4), in order to Increase the mechanical 
strength of the assembly and to decrease windage losses. This assembly 
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hes been spin tested by Delco at speeds greater than 30,000 rpm without 
failure. This is greater than three times the rated speed of 9,000 rpm. 
The small size of this machine is readily apparent in the actual sized 
cross section of the motor shotro in Figure (2.1-5). 

The stator assembly consists of a stack of 7 mil Vanadium Permendur 
punchings, see Figure (2.1-6). Each stator punching consists of 24 slots 
and 24 teeth. Figure (2.1-7). The slots are only partially filled in 
order to allow the cir<"ilation of the coolant. A nonmagnetic sleeve 
covers the inside of the stator bore thereby isolating the stator cooling 
system as well as reducing windage losses. 

A mechanical brake is provided with each motor in order to prevent 
rotor motion whenever the machine is turned off. This brake makes use 
of a solenoid actuated braking mechanism, see Reference 

At the other end of the rotor shaft lies the rotor position sensor 
or RPS which is used to control the switching of the inverter transistors. 
The RPS consists of a laminated rotor, containing four lobes, which 
operates in conjunction with an electrical sensor composed of nine coils 
on a laminated stator. The nine coils are arranged into three groups of 
three coils each. The resolution of the RPS is 60 electrical or 15 
mechanical degrees. The logical output signals of these three coil groups 
(i.e. AA, BB and CC) versus rotor angle are shown in Figure (2.1-8). 

The RPS output at some angle 0 can be simulated by comparing to zero the 
values of the three sinusoids shown as dashed lines in Figure (2.1-8) 
as follows: 


^ fO IF SIN[4(0g + 0)3 < 0 

IF SIN[4(©s + 0)] 1 0 


( 2 . 1 - 1 ) 
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and 



IF SIN[4(0g +0 - 30*)] < 0 

IF SIN[4(0g + 0 - 30*)] ^ 0 


0 IF SIN[4(0g + 0 - 60*)] < 0 

1 IF SIN[4(0g + 0 - 60*)] ^ 0 


( 2 . 1 - 2 ) 


(2.1-3) 


The angle 0 represents a constant shift In the output of the RPS. Figure 
s 

(2.1-8) shows the RPS output versus rotor angle for 0 ■ 0*. 

s 

A tachometer is attached to that end of the rotor shaft which is 
opposite to the end containing the RPS. The output voltage of the tach- 
ometer versus motor speed follows a gradient of 7.0 volts per 1000 revo- 
lutions per minute. This output signal is fed back to the channel control 
electronics where it is processed along with other data to generate the 
EMA control signals. 
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2.1.2 Electrical Description 

The permanent magnet brushless d.c. motor Is rated at 270 volts 
peak line to line and 60 amperes phase current at rated speed and torque. 
The electrical configuration of this machine is determined by the eight 
pole permanent magnet rotor assembly described in the previous section. 
Consequently, the relationship between electrical and mechanical angles 
is given by: 

- P© /2 
e m 

or 

©^ - 4© (2.1-4) 

e m 

where 

P • Number of Poles 

©g » Electrical Angle 

© * Mechanical Angle, 

m 

The stator has 24 slots which means that there are six slots per 
pole pair or three slots per pole pitch and two slots per phase. Figure 
(2.1-9) gives a developed view of a cross section of the stator and rotor. 
The angular reference for the stator is taken from the phase (a) MMF axis 
farthest from the right hand side. The rotor angular references are 
taken at the centers of the north poles. 

The structure of the stator phase windings per pole pair is shown 
in Figure (2.1-10). Each phase winding consists of a series connection 
of four sets of two-coil groups connected in parallel. Each coil group 
consists of 11 turns of two AWG #18 heavy ML wires connected in parallel 
which results in a total of 44 AWG #18 HML conductors per slot. However, 
there are only 11 independent conductors per slot since the 44 conductors 
are connected into 11 groups of four conductors in parallel. These 11 





groups of conductors are equivalent to 11 turns per pole pair of a con- 
ductor with a cross sectional area equivalent to four AWG //18 HML con- 
ductors. The phase windings of the four pole pairs are connected together 
In series resulting in an equivalent phase winding of 44 turns. 

The open circuit voltages at various speeds, as determined through 
tests performed by Delco, Reference (2), are given in Table (2.1-1). 

A typical open circuit voltage waveform is shown in Figure (2.1-11). 

Based upon Table (2.1-1) the relationship between the open circuit voltage 
between the open circuit voltage and the rotor speed is given approximately 
by: 

E - 0.01555 n (2.1-5) 

oc 

where 

E ■ Peak open circuit phase voltage in volts 

OC 

n ® rotor speed in rpm 

The electrical operating modes or regimes of this machine are 
analogous to those of conventional synchronous machines. These operating 
modes are: motoring, regenerative braking, and plugging. The configura- 

tions of this machine and the power conditioner during these modes are 
the topics of the next section. 
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2.2 MACHINE MODES OF OPERATION 

The EMA motors during an actuation maneuver can function In any 
one of three possible operating modes. The operating regions of these 
modes In the torque-speed plane are shown In Figure (2.2-1). The modes 
outlined In this figure are: 

a) The Motoring Mode - Quadrants I and III, 

b) The Regenerative Braking Mode - Quadrants II and IV, and 

c) The Braking Mode - Quadrants II and IV. 

An analysis of these operating modes Is the topic of the three sections 
which follow. This analysis Is based upon the sample data presented in 
Table (2.2-1) which was obtained from the VPI EMA computer model. 
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Table (2.2-1) The Status of the Inverter during the Motoring, Regenera- 
tive Braking and Plugging Modes at a Rotor Angle of 16.7®. 


I. Motoring Mode 

Speed • 4000 rpm Rotor Angle ■ 16.788® (mech) 

01 - ON, 02 - OFF, Q3 - OIF, Q4 - OFF, Q5 • OFF, Q6 - ON 
D1 - OFF, D2 - ON, D3 - OFF, D4 - OFF, D5 • OFF, D6 - OFF 
Phase A Current ■ 60.67472841 Amperes 
Phase B Current • -23.76844021 Amperes 
Phase C Current ■ -36.90628439 Amperes 

II. Regenerative Braking Mode 

Speed ■ 4000 rpm Rotor Angle ■ 16.788® (mech) 

Q1 • OFF, Q2 » OFF, Q3 ■ OFF, Q4 ■ OFF, Q5 • OFF, Q6 - OFF 
D1 - ON, D2 « OFF, D3 » ON, D4 • OFF, D5 - ON, D6 » OFF 
Phase A Current * -22.34997178 Amperes 
Phase B Current » 25.27216082 Amperes 
Phase C Current ■ -2.922189036 Amperes 

III. Plugging Mode 

Speed - 240 RPM Rotor Angle « 16.740® (mech) 

Q1 » OFF, Q2 » ON, Q3 * OFF, Q4 « ON, Q5 » OFF, Q6 * OFF 
D1 » OFF, D2 = OFF, D3 - ON, D4 - OFF, D5 » ON, D6 OFF 
Phase A Current » -44.47176035 
Phase B Current « 62.56684007 
Phase C Current * -18.09507972 

*This data was obtained from typical computer simulations. 
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2.2.1 Motoring Mode 

The motoring mode Is used to accelerate the system inertia in the 
direction of motion, or in other words, the motor output torque is in the 
direction of motion (quadrants I and III as shown in Figure (2.2-1)). 

There is, of course, a practical limit to the amount of torque and speed 
which can be produced in this mode because of the following three factors: 

a) The chosen power supply voltage limits the maximum speed at which 
the machine can deliver the rated torque. This is due to the in- 
crease in the phase back emfs with increases in speed. These back 
emfs oppose the power supply voltage. Consequently when the 
speed is such that the two voltages are equal, the current can 

no longer be driven into the machine phase windings resulting 
in reduced torque and speed. 

b) The volt-ampere rating of the switches limits the amount of power 
which may be Injected into the machine phase windings. 

c> The volt-ampere rating of the machine itself sets the upper bound 
on machine torque. 

The maximum obtainable torque in this machine is basically limited 
by the volt-ampere rating of the switches (item (b) above). This maxi- 
mum torque is symbolized by the constant TMMAX in Figure (2.2-1). Pre- 
sently, TMMAX is set at 9.0387 nt-m (80 in- lb) which is approximately 
equivalent to a phase current of 40 amperes. 

In order to gain Insight into the operation of this electromechanical 
act\iator one must examine the behavior of the power conditioner switches 
as well as the interaction of the electromagnetic fields inside the 
motor. Figures (2.2-2) and (2.2-3) display typical current paths in 
the power conditioner with the chopper transistor QM on and off respec- 
tively and with Q1 and Q6 on, see item (I) In Table (2.2-1). The chopper 
transistor QM switches on and off in order to maintain the current in 
the coupling inductor within a specified band about a commanded current 
level. With QM on. Figure (2.2-2), the battery is connected to the 
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coupling Inductor which results in an increase in the inductor current. 
Transistor QM is turned off when this current reaches a specified level 
above the commanded current > see Figure (2.2-3). The current then de- 
creases until the inductor current drops to a specified point below the 
commanded current; at which point QM is switched on again thereby repeating 
the cycle. 

The controlled coupling inductor current is then injected into the 
proper phase windings of the machine via the six transistor and six diode 
inverter bridge (Q1 - Q6 and D1 - D6). The firing of these transistors 
is controlled by the output of the rotor position sensor or RPS. The 
logic which translates the RPS signals into commands to switch the in- 
verter transistors on or off is described in detail in Section (4.2.4) 
and in Reference (2). 

As an example, consider the two typical network configurations 
shown in Figures (2.2-2) and (2.2-3). In both cases, Q1 and Q6 are on 
which means that coupling inductor current is flowing into phase (a) 
and out of phase (c). These two phase currents result in a motor torque 
which accelerates the system inertia in the direction of motion. The 
current out of phase (b) is in the process of decaying through diode D2 
because of the recent cutoff of transistor Q5. Diode DR is also on 
in this case in order to provide a path for the excess current resulting 
from this commutation process. 

The electromagnetic interactions inside the permanent magnet orush- 
less d.c. motor for the case shown in Figure (2.2-2) and listed under 
item (I) of Table (2.2-1) are displayed in Figure (2.2-4) for a single 
turn of phase (a) . The symbols used in this diagram are defined as 


follows : 







VBAT ■ Battery voltage (volts) 
f m Conductor Current Vector (Amperes) 

* o 

* Conductor Flux Density Vector (Webers/Meter ) 

V^ « Velocity Vector of the Conductors with respect 
to (Meters /Second) 

E • Induced Conductor Back Emf Vector (Volts /Meter) 

■ Conductor Force Vector (Newtons/Meter) 

Fj^ * Rotor Force Vector (Newtons) 

Vj^ • Rotor Velocity Vector (Meters/Second) 

The force per unit length, F^, which is experienced by a conductor having 
the current vector, I, and a flux density vector, B^, Is given by 

F^ - (I X 5^)1^ (Newtons) (2.2-1) 

where Is the active conductor length. For the two conductors shown In 
Figure (2.2-4) the forces F^, are both In the negative direction. Since 
the stator Is fixed and according to the action and reaction principle, 
the sum of these conductor forces Is translated Into a rotor force, F^, 
which Is In the direction of rotation and Is opposite In direction to the 
rotor conductor force vector, F^. The orientation and magnitude of the 
Induced back emf, E, per unit length for a conductor experiencing a flux 
density B^(wb/m ) and having a velocity V^ (m/s) with respect to the field 
Is given by the vector product 

E » V- X B-. (Volts /Meter) (2.2-2) 

c c 

In the example shown In Figure (2.2-4), the Induced back emf opposes the 
flow of the phase current, I, as would be expected In the motoring mode. 
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2.2.2 Regenerative Braking Mode 

The regenerative braking mode utilizes the fact that these per- 
manent magnet machines when operated without the power conditioners look 
like salient pole synchronous generators. During this mode, the permanent 
magnet machine Is operated as a generator thereby converting the energy 
of the rotating masses Into electrical energy. This energy Is used to 
control the braking torque as well as to recharge the supply batteries. 

In this mode the rotor torque and velocity vectors point In opposite di- 
rections which results In a deceleration of the rotating masses. This 
relationship between torque and speed Is graphically displayed In Figure 
(2.2-1). The regenerative braking mode Is confined by the control sys- 
tem to the speed regions designated by SFDLO and SPDHI In Figure (2.2-1). 
This Is necessary because outside of these speed regions the machine Is 
turning too slowly to produce the back emf necessary to maintain the com- 
manded braking torque (line current). 

When the actuator motors are operated In this fashion all six In- 
verter switches (Q1 - Q6) are turned off. Consequently, the phase 
currents must pass through the six inverter diodes (D1 - D6) which to- 
gether act as a full wave rectifier bridge. Figures (2.2-5) and (2.2-6) 
show some typical current paths during this mode of operation, with QB 
on and off respectively (see item (II) in Table (2.2-1). In Figure (2.2-5) 
diodes Dl, D3 and D5 are on and the chopper transistor, QB, Is off. Con- 
sequently, the line and phase currents will increase. With QB on, the 
battery receives no current from the machine. This changes as soon as 
QB is turned off. When QB is off. Figure (2.2-6), current flows from the 
inverter bridge through the braking diode, DB, imd into the battery. In 
this configuration the EMA actually converts the externally applied load 
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torque into electrical energy which is used to recharge the system 
batteries. 

The electromagnetic Interactions between a single turn of a stator 
phase winding carrying a current, I, and one north-south pole pair are 
shown in Figure (2.2-7) for the case given in Figure (2.2-6). The rotor 
velocity V^, the Induced back emf E, and the flux density, B^, all have 
the same orientation as in the motoring mode, Figure (2.2-4). The dif- 
ference lies In the direction of the phase current I. In Figure (2.2.7) 
the phase current exits through the positive terminal of phase (a) and 
enters the positive terminal of the battery. In other words, the stator 
winding acts as a source of power while the battery acts like a load. 

This Is opposite to the motor case discussed earlier. The majority of 
the energy received by the battery is stored and can later be made availa- 
ble. The negative phase current results in a conductor force, F^, which 
points in a direction opposite to V^. Consequently, the rotor force 
vector, Fj^, and the rotor velocity vector, point in opposite di- 

rections which results in a deceleration of the system inertia as expected. 
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2.2.3 Plugging Mode 

The EMA operates in the plugging mode whenever the motor speed is 
too low for the regenerative mode to supply the commanded braking torque. 
The operating regions for the plugging mode in the torque-speed plane are 
shown in Figure (2.2-1). 

Figures (2.2-8) and (2.2-9) show some typical current patterns in 
the power electronics during the plugging mode with the chopper transis- 
tor QM on and off respectively. Notice that in both cases Q2 and D5 are 
both on. This effectively sets the voltage of the node formed by the 
Intersection of the coupling inductor and the switches DR, Ql, Q2, Q3, 

Dl, D2 and D3 close to ground potential. This switch configuration 
allows circulating currents in the inverter switches and machine phase 
windings, see loops and I 2 Figure (2.2-8). These currents do not 
pass through the chopper transistor or coupling inductor and consequently 
they cannot be controlled by the EMA control system. For this reason 
the current in the coupling Inductor is limited to 25 amperes in order to 
keep the currents in the inverter switches below the 60 ampere maximum. 

This limit on in the maximum chopper current results in a reduced maxi- 
mum braking torque which is symbolized by TPMAX in Figure (2.2-1). 

The electromagnetic interactions between the rotor and stator cir- 
cuits during the plugging mode are depicted in Figure (2.2-10) for the 
case listed under item (III) iiv Table (2.2-1). The plugging mode looks 
almost identical to the regenerative braking mode. Figure (2.2-7), in terms 
of the electromagnetic interactions inside the machine. The difference 
lies in the polarity of the battery with respect to the induced back emfs. 
In the regenerative braking case. Figure (2.2-7), the battery is connected 
so that it opposes the machine back emfs. Consequently, rotational energy 
of the rotating mass is converved into electrical energy. In the plugging 
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mode, on the other hand, the battery voltage aids the machine back .tmf.?. 
Figure (2.2-10), whicii means that the battery supplies energy to sl» w 
down the rotating masses. 


» 
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2.3 ELECTRICAL MODELS OF THE PERMANENT MAGNET MACHINE 

Two different electrical models of the permanent magnet brushless 
d.c. motor have been considered as possible candidates for the computer 
model developed under this contract. The first model consists of the 
standard six winding machine model commonly used to analyze conventional 
synchronous machines while the second model Is based upon a back emf be- 
hind a leakage Inductance. The back emf behind a leakage Inductance 
machine model was chosen because of computational considerations and 
because the parameters of such a model (l.e. back emfs and leadage In- 
ductances) can be obtained readily from a nonlinear finite element field 
solution. These models are described in greater detail in Sections 
2.3.1 and 2.3.2 which follow. 
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2.3.1 Six Winding Machine Model 

As stated earlier, the structure of these permanent magnet brushless 
d.c. motors Is very similar to the structure of conventional salient pole 
synchronous machines. Consequently, one can apply existing state space 
models of conventional synchronous machines to this new class of permanent 
magnet machines. The only difference between these machines and conven- 
tional salient pole synchronous machines Is the lack of the field and 
damper windings on the rotor. The field winding function Is performed by 
the permanent magnet pole pieces. 

The six winding machine model developed In this section Is based 
upon the lumped circuit model of a conventional two pole salient-pole 
type synchronous machine as shown In Figure (2.3-1). This lumped circuit 
model consists of three stator phase windings (a, b, c), two damper wind- 
ings (kd, kq), and a field winding (f ) . A schematic representation of 
these windings Is given In Figure (2.3-2). 

The field winding shown In these figures does not actually exist In 
the machines at hand since permanent magnets are used to provide the 
necessary magnetic field. However, a "fictitious" field winding Is re- 
tained In this model as an artifice to represent the magnetomotive force 
produced by these magnets. In order to obtain a value for this "fictitious" 
field current, ip, one uses the demagnetization curve of the permanent 
magnet. Figure (2.3-3). Multiplying the coercivity, (Amperes /Meter ) , 
by the height of a magnet, (Meters), one obtains the field mmf, 

F^^ , . ■ H *1 (Ampere-Turns) (2.3-1) 

field c m 

This is equivalent to a fixed value of field current for a given number 


of fictitious turns. 
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The six winding machine model can now be derived by writing the 
differential equations corresponding to these windings. Consider, for 
example, winding (a) in Figure (2.3-2) which is the stator phase winding 
corresponding to phase (a) . Assuming that this winding has a resistance 
of R (Ohms) the phase voltage, v , can be written as 


R i + 77 X 
a a dt a 


where 


and 


''a ■ *a^a + It * ^abS * ^ac^c 


+ + L. 


“af*f ' “akdStd ^akq^q 


(2.3-2) 


X ■ Total flux linkage of coil (a) (Weber-Turns /Ampere) 

L ■ Self inductance of coll (a) (Henries) 
aa 


^ab* ^ac* ^af* ^akd ^akq mutuals between coil (a) 


'' and coils b, c, f, kd and kq respectively. (Henries) 

For phases (b) and (c) similar expressions are written as follows: 


and 


dt ^b 


\ ‘ dt ^^ab^a SbS Hc^c 

■** h>f^f Skd^d ■*' HkqScq^ 


R i + X 
c c dt c 


(2.3-3) 


v-Ri +77 (L i + L.i,+L i 
c c c dt ac a oc d cc c 


^cf^f ^ckd^d * ^ckq^q^ 


(2.3-4) 


S imi lar expressions are written for the coils of the field and damper 
windings f, kd, and kq respectively. The resulting system of six coupled 
first order differential equations constitutes the six winding lumped 
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circuit parameter model of a synchronous machine. In matrix form these 


equations are written as follows: 


r 1 

V 

a 


a)’** 

lb 

o 

o 


"^a 1 



Ol 

o 

1 

o 



V 


0 0 R 1 


i 

• £ - 


c 


. £ - 

^f 


1 Rj 0 0 


if 

0 


O 

0 

01 


ikd 

0 


1 

cr 

o 

o 

1 


Aq. 


i'aa ^ab ^ac [ i*af i*akd ^akq 


ia 

N 

i*ab Sb h)c • Sf Hkd Skq 


ib 


i'ac ^hc ^cc I ^cf i'ckd ^ckq 




i*af Sf ^cf 1 ^ff ^fkd ^fkq 




i*akd Skd ^ckd ' ^fkd \dkd ^dkq 


^kdl 


^akq S)kq i*ckq , ^fkd \dkq \qkq 

m 


Sc,: 

m fli 

/ 


All of the inductance coefficients in equation (2.3-5) above are de- 
pendent upon the angular position of the rotor and upon the magnetic satu- 
ration picture which is a function of the currents. Hence all the induc- 
tance coefficients are functions of time. Therefore the inductance matrix 
cannot be taken outside the derivative operator (-^) . 

The samarium-cobalt permanent magnet pole pieces offer a high resis- 
tance to the flow of eddy currents because of their low electrical con- 
ductivity (0.2 X 10'' o/m versus 5.8 x lo^ O/m for copper). Consequently 
the currents i^^ and 1^^ which flow in the fictitious (d) and (q) axis 
damper windings kd and kq are negligible. Applying these simplifications 
to Equation (2.3-5) reduces the order of this machine model from sixth to 


fourth order as follows: 
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+ 


dt 


aa 

^b 

1 

^af 

ab 

Hb 

He : 

Sf 

ac 

mm 

^Hc_ 

He i 

^cf 

■af 


h;i 

Lff 


a 

s 

(2.3-6) 


Equation (2.3-6) represents a lumped parameter model of the per- 
manent magnet brushless d.c. motor at hard. Unfortunately this model is 
rather combersome because of the fact that the inductances are all func- 
tions of time. There are techniques for determining the values of these 
inductances. Reference (2.3-1). However, a simpler but quite adequate 
approach to modeling this machine was adopted and is the topic of the 
next section. This adopted approach has proven, itself, after comparing 
its numerical results with experimental test data, to be very accurate 
and reliable in the prediction of the dynamic performances of the machine. 
Proof of this is given in later chapters of this report. 
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2.3.2 Transient Back Emf - Leakage Inductance (Transient Inductance) 

Machine Model 

The back emf-leakage Inductance machine model per phase consists of 
a series connection of three elements; the d.c. winding resistance, the 
winding leakage (transient) Inductance, and the Induced winding (transient) 
back emf. A schematic diagram of this model Is shown In Figure (2.3*4) 
for the three phase machine under study. This model will be derived from 
the reduced six iflndlng machine model of equation (2.3-6) In the develop- 
ment which follows. 

The stator phase inductances and of equation (2.3-6) 

can be expressed as the sum of two Inductances: 


^aa “ 

^la *'■ ^ga 

(2.3-7) 

XT 

II 

^Ib * ^gb 

(2.3-8) 

II 

U 

o 

^Ic ■*’ ^gc 

(2.3-9) 


That Is, a leakage Inductance associated with the leakage flux and an In- 
ductance associated with that portion of the total flux which crosses the 
alrgap, see Figure (2.3-5). These inductances are symbolized by 
^Ib, ^Ic^ ^^a» ^b» ^c^ respectively and they are defined as fol- 
lows: 


la 


'lb 


Ic 


ga 


gb 


Vb'S 

Ic c 


“♦itb/S, 

N* /i 
ga a 


X /I • N, /I 

gc C gc C 


L 

gc 



gc 




ge and Air Gap Flux Paths 







where 


X , X., , and X are the leakage flux linkages of phases 
® (a) , (b) and 'c) respectively; 

X , X , and X are the flux linkages of phases (a), (b) 

8® 8^ 8C and (c) respectively, corresponding to 

that portion of the flux which crosses the 
air gap; 

(|>- , , and are the leakage fluxes of phases (a), (b) 

(c) respectively; 

^ea* ^eb» ^ec fluxes of phases (a), (b) 

and (c) respectively; 


N equals the number of turns of the abc phase windings. 

Substituting equations (2.3-7), (2.3-8) and (2.3-9) into equation 
(2.3-6) results in: 

■ R 0 0 

a 

0 R^ 0 

0 0 R 

_c 

0 0 0 


^ i 

dt ' 
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(2.3-10) 
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or in symbolic matrix notation one can rewrite equation (2.3-10) as follows 


•^bcf ^bcf-^bcf ^ dt ^^bc£ ^bcf 


dt ^^bcfg ^bcf^ 


(2.3-11) 


The slot leakage shovm in Figure (2.3-5) is mainly confined to an area 
immediately surrounding the slot conductor, therefore much of the leakage 
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flux passes through nonmagnetic materials such as copper , air , and Insu- 


lation which do not experience saturation. In addition, the operating 

2 2 

air gap flux densities at hand, namely, 0.64 wb/m (41000 Llnes/ln ) are 
relatively low, which reduces the levels of saturation In this machine. 
Consequently, the leakage Inductance depends mainly upon the slot geometry 
and to a lesser degree upon the state of the machine currents. For these 
reasons, the leakage Inductance matrix, essentially time Invarl- 

/ d V 

ant and can therefore be moved outside of the derivative operator, » 
equation (2.3-11) then becomes: 


^bcf " ^bcf-^bcf ^ ^abc£ dt ^-^bcf^ 
^ ^^bcfg-^bcf ^ 


(2.3-12) 


The calculation of the "ficticious” field winding voltage, v^, in 
equation (2.3-10) is unnecessary in the case of the permanent magnet ma- 
chine at hand because it has no physical existance. Consequently, the 
fourth row of equation (2.3-10) caa be eliminated as follows: 
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(2.3-13) 
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or In compact notation 


^bc ^bc ^bc ^ ^bcl dt ^bc 


■*' d^t ^^%bcfg -iabcf^’ 


(2.3-14) 


The three right hand terms of equation (2.3-14) are nothing more 
than the components of the back emf - leakage inductance machine model 
shown earlier in Figure (2.3-4), where 


~abc 


R 0 
a 


(2.3-15) 


represents the stator phase winding resistances, 


^bc£ 


ha ° ° 

° hb 0 

0 0 L. 


(2.3-16) 


represents the stator leakage inductances, and 
■^bc " dt ^hbefg ■^bef^ 


(2.3-17) 


represents the phase back emfs, where e^^^ is the vector containing the 
phase back emfs; e^, e^, and e^. 

Examination of equations (2.3-14) and (2.3-17) reveals that for the 
case where i^j^^ ■ 0_, the terminal voltage vector, becomes: 


V = e 
— abc -abc 


^af 

hf h 


(2.3-17) 


In other words the no load or open circuit teminal (phase to neutral) 
voltages are equivalent to the no load or open circuit back emfs. The 
open circuit back emfs can be obtained either through an open circuit 
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test as Is the case In Figure (2.1-11) or by means of a finite element 
solution. Regardless of how the open circuit emf waveform is determined, 
it must be modified whenever stator currents are present. This can be 
accomplished by calculating the resultant air gap MMF distribution due to 
the nonzero ■ stator currents and then magnetizing or demagnetizing the 
magnetic pole pieces of the rotor accordingly. These back emf calcula- 
tions are the topic of the next section. 
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2.4 CALCULATION OF MACHINE BACK EMFS AND TORQUE 

The machine model derived In the previous section requires the 
knowledge of the air gap phase back emfs on an Instantaneous basis. 

The back emfs are functions of: (a) rotor position, (b) rotor speed, 

(c) the stator phase currents, and (d) the magnet properties. 

A procedure to approximate the back emfs as a function of these 
parameters Is presented In this section. In addition, the Instantaneous 
machine torque Is calculated from the values of the Instantaneous phase 
back emfs, phase currents, and machine velocity. The calculated machine 
torque later becomes an Input (or forcing function) to the dynamic model 
of the rotating masses which is developed In Chapter 4 of this volume. 
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2.4.1 Calculation of Machine Back Emfs 

The determination of the machine back emfs requires detailed 
knowledge of the flux density distribution In the air gap at any given 

® Hot* 

Instant of time. Such distributions are obtainable from an Instantan- 
eous finite element solution of the nonlinear Poisson's equation which 
governs the machine behavior. References (7) and (8). The problem with 
this approach Is the computational effort required to solve the large 
systems of nonlinear algebraic equations which are produced by the 
finite element discretization at each Instant of time. For this reason, 
a lumped parameter model of the machine magnetic circuit was adopted. 

The lumped parameter approach approximates the distributed nature of 
the magnetic fields by a lumped parameter network consisting of equiva- 
lent flux path reluctances and magnetomotive forces or mmfs. 

For purposes of discussion, a developed cross sectional view of 
the magnetic circuit spanned by two pole pairs (180® mechanical, 720® 
electrical) is given in Figure (2.4-1). This figure displays the 
electromagnetic interactions between two pole pairs of the stator and 
rotor circuits under maximum torque conditions. Maximum machine torque 
is developed when the stator and rotor mmfs are displayed by 90 elec- 
trical or 22.5 mechanical degrees. References (13) and (14). The 
stator phase winding which produces the stator mmf is represented by 
equivalent current filaments located at the centroids of the stator 
slots. These conductors are identified by the letters A, B, and C 
corresponding to phases, and by a superscript of (+) or (-) depending upon 
the orientation (+ or - z) of the induced back emf, for a positive 
phase current during motoring in the positive direction. 
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The simfs FA, FB, and FC produced by the above mentioned phase 
currents are displayed below the developed machine cross-section, Figure 
(2.4-1). The resultant stator mm£, FS, Is obtained by summing these 
phase mmfs. The Idealized mmf waveform, FM, of the permanent magnet 
pole pieces Is given directly below the phase mmfs. These mmfs estab- 
lish magnetic flux which follows closed paths similar to the ones shown 
In the machine cross-section. 

The electromechanical Interactions between the stator and rotor 
can be approximated by a lumped magnetic circuit consisting of reluc- 
tances In series with sources of mmf. Applying this model to any one 
of the eight symmetrical closed flux paths between the stator and rotor 
circuits yields the magnetic network of Figure (2.4-2). The permeability 
of the iron rotor and stator cores can be assumed to be Infinite because 
of the dominance of the air gap and magnet reluctances at the low levels 
of flux densities In these cores. This model of the machine forms the 
basis from which the phase back emfs are calculated. 

As was shown earlier in Section 2.2, one can determine the in- 
duced emf in a phase conductor by means of equation (2.2-2), once the 
flux density, B^, at the conductor and the velocity, V^, of the con- 
ductor with respect to the field, B^,, are known. In the air gap, B^ 

Is essentially radially oriented, and therefore perpendicular to V^. 
Consequently, the magnitude of the induced emf per unit conductor 
length, E, as defined by equation (2.2-2) becomes 

1E| - |V^| • |Bj (2.4-1) 

Before the back emfs can be calculated, the total number of 
effective series conductors per phase must be determined. Each phase 
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winding per pole pair consists of 11 effective turns. Since there are 
four pole pairs, the total number of turns per phase winding equals 44. 
Each of these turns is composed of two conductors which experience equal 
but oppositely oriented flux densities, see Figure (2.2*'4). Therefore, 
the magnitude of the emf per unit length, |e| , which is induced in each 
of the 88 conductors is equal. The total phase back emfs Induced in the 
phase windings therefore become: 






Volts 

(2.4-2) 

®B ■ 




Volts 

(2.4-3) 

®C ‘ 

188 VzI • 


- 'c®sc» 

Volts 

(2.4-4) 


All parameters in equations (2.4-2), (2.4-3), and (2.4-4) are 
known, except for the flux densities 

These flux densities are functions of the rotor position, stator cur- 
rents, and of the magnet demagnetization characteristics. In order to 
keep the computational costs within manageable bounds, the just mentioned 
lumped parameter approach will be used. In this approach, the flux 
path around the phase conductors is replaced by an equivalent lumped 
parameter circuit as shown in Figure (2.4-2). This circuit consists 
of an iron core, representing the flux paths in the stator and rotor 
bodies connected in series with an air gap and permanent magnet. 

Notice that only one half of the total flux path shown in the developed 
diagram of Figure (2.4-1) needs to be considered because of symmetry. 

The behavior of the magnet is modeled by the magnet demagneti- 
zation (B-H) curve. The B-H characteristic of the permanent magnets 
used for the rotor pole pieces is given in Figure (2.4-3). The upper 
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line represents the INTRINSIC demagnetization curve of the permanent 
magnet and in the linear region is described by the equation: 

J - -(l*„ + (2.4-5) 

The lower line is called the NORMAL demagnetization curve and in the 
linear region can be approximated by the line: 

‘m • + h 

These two curves obey a linear interrelationship given by 

. J + (2.4-7) 

As mentioned before, the flux paths in the stator and rotor iron are 
assumed to be ideal (y ■ «) because of the relatively low operating 
flux densities encountered in this machine. 

At the present time, the air gap plus magnet length, is 
held fixed for all levels of machine torque. However, a finite ele- 
ment analysis would allow the calculation of an effective air gap 
length based upon different levels of machine torque. As the torque 
increases, the flux is twisted, thus increasing the effective air gap 
length, see Figure (2.4-1). 

The procedure used to calculate the flux densities of the stator 
conductors is as follows: 

Step 1 : Determine the open circuit (stator phase currents zero) 

/ 

flux densities from a Fourier analysis of the open circuit voltage 
(back emf) waveform given in Fig ire (2.1-11). A normalized series 
representation of this waveform is given by: 


J 
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1.111111 cos Sj^ 

- .0819017 cos 3Gj^ 

- .0807241 cos (2.4-8) 

+ .0405518 cos 

This series Is multiplied by the flux density, BPEAK, such that substi- 
tution Into equation (2.4-2), (2.4-3), and (2.4-4) results in the mea- 
sured open circuit voltage waveform given In Figure (2.1-11). The above 
series therefore becomes: 


BGAP » BPEAK [1.111111 cos 

- .0819017 cos 30jjj^ 

- .0807241 cos 50j^j^ (2.4-9) 

+ .0405518 cos 70j^ 

The spatial flux density distribution represented by this series Is 
produced entirely by the rotor permanent magnets. Rewriting this series 
in terms of a stator reference angle 0^^ and the angle 0^^ between the 
phase A mmf axis and the closest rotor north pole, Figare (2.1-9), yields 


BGAP - BPEAK [1.111111 cos [4(0gg - Ogj^)] 

- .0819017 cos 3 [4(0-5 - ©- ) ] 

- .0807241 cos 5[4(©gg - ©gj^) ] (2.4-10) 

+ .0405518 cos 7[4(©gg - ©gj^) ] 

f 

This equation ' enables one to determine the air gap flux density at any 
stator angle 0gg given a rotor-stator displacement of ©g^^. The open 
circuit conductor flux densities can novi be calculated simply by sub- 
stituting for 0gg the angular displacements ©g^, ©gg, and ©g^ of the 
conductors. 


Step 2 . Determine the equivalent magnet and air gap lengths 


and £g necessary to produce the flux densities given in STEP 1. 
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These calculations make use of the magnetic lumped circuit machine 
model of Figure (2.4-2) and the magnet demagnetization curve of 
Figure (2.4-3). In the open circuit case the stator nanf, FS, equals 
zero, therefore by Ampere's Law one obtains: 

/H-dl-HA +H1 -0 (2.4-11) 

m m g g 

Since fringing Is neglected 

B “B *B*uH (2.4-12) 

mg 0 g 

Substituting equation (2.4-12) into (2.4-11) and solving for results 
In: 

H • (2.4-13) 

ID g 0 m 

Substituting equation (2.4-13) into equation (2.4-6) yields 

B = (BT,l„/y„£ Hp)B + B^ (2.4-14) 

K g 0 HI c r 

The ratio i H can now be written as 
g m 

i.Jl = (2.4-15) 

g ID O C K O L 

Now 

1 +1 (2.4-16) 

1 g ni 

Substituting for in equation (2.4-15) and solving for results 
in: 

S ■ - B^] + Bj^B 


(2.4-17) 
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Step 3. Determine the resultant stator mmf, Fg, at conductor 
locations. The stator mmf at any arbitrary stator angle, 6gg, can be 
approximated as shown in Figure (2.4-1). 


Step 4. Given 2^ and Fg determine the new operating flux 
density B'. From H • d2 ^ I one obtains: 


Solving for H results in; 

HI 


H 


ffl 


H„2m + H^2 - F- 

mm g g s 


7JI - B'2^/u 2^ 
5 m g o m 


(2.4-18) 


(2.4-19) 


Substituting equation (2.4-19) into equation (2.4-5) yields the intrin- 
sic flux density J as follows; 


J - \ <2.4-20) 

The new operating flux density B' of the magnet is given by Equation 
(2.4-7). Substituting for and J into equation (2.4-7) and solving 
for B' yields; 


B' = 


- F3I 

“c'm - 


(2.4-21) 


This flux density is then used to calculate the induced back emf as 
described earlier. In the subsection which follows, these emfs are 
used to determine the developed electromagnetic machine torque. 


2 . 4.2 Machine Torque Calculation 

The developed machine torque is required as an input for the 
model of the dynamics of the rotating masses which is developed later 
in Chapter 4. The developed air gap electromagnetic torque is calcu- 


lated as follows: 
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TEM • (e^i^ + ejij + fej,i^)/RVEL (2.4-22) 

where TEM Is the developed air gap electromechanical torque; 
e^, Cg, and are the phase back em£s; 

1^, Ig, and are the phase currents; 

and RVEL is the rotor velocity. 


i 
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3. STATE SPACE MODEL OF THE POWER CONDITIONER AND MACHINE NETWORK 

The utilization of the permanent magnet machine, described in the 
pervious chapter, for servo applications requires a current or power 
conditioner network capable of controlling the machine phase currents 
so as to produce actuation torques of specified orientations and 
magnitudes. This function is performed by four power conditioner 
networks, one for each of the four machines ;fhlch together constitute 
a single four channel EM actuator. 

The nonlinear differential equations and their related output 
equations which characterize the electromagnetic dynamics of a single 
power conditioner and machine unit are derived in this chapter. The 
differential equations are •subsequently Integrated forward in time 
using the exponential series method. This projection of the system 
state variables forward in time makes it possible to calculate the 
developed machine torques as well as the power conditioner voltages, 
currents, and powers instantaneously. This time sequence of machine 
torques couples the electromagnetic model of the machine and power 
conditioner developed in this chapter with the mechanical model 
developed in the next chapter. The calculated, instantaneous 
values of voltages, currents, and powers are used to quantify the 
stresses experienced by the power conditioner components under various 
operating regimes. The knowledge of the magnitude and nature of these 
stresses is of great concern in aerospace applications where system 
reliability and margins of safety are of paramount importance. 
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3.1 DESCRIPTION OF THE POWER CONDITIONER 

The function of the power conditioner as stated earlier in the 
introduction to this chapter is to modulate the d.c. current of the 
supply battery into a three phase current which will produce a specified 
average machine torque in response to a torque (current) command from 
the servo controller. Figure (3.1<-1) shows a schematic diagram of 
the power conditioner along with the associated three phase permanent 
magnet machine. The principal components or elements of the power con- 
ditioner can be categorized as follows: 

a) A DC Power Supply 

b) A Filter Capacitor 

c) A Pulse Width Modulator or Chopper 

d) A Coupling or Chopper Inductor 

e) A Six Leg Power Inverter Bridge 

These power conditioner elements are indicated in Figure (3.1-1). 

The d.c power supply, in the present system, consists of a bank of 
24 conventional lead acid batteries ccunected together to produce a 
nominal potential of 270 volts. This is roughly equal to the d.c. voltages 

which presently exist on the space shuttle. The battery output is 

connect ''d across a large filter capacitor which helps to suppress 
voltage transients which might otherwise result in damage to the switch- 
ing transistors QM and QB as well as suppressing possible supply 
transients. 

The d.c. power available at the terminals of this capacitor is then 
pulse width modulated or chopped into a series of current pulses. The 

width of these pulses is controlled by logic described in the next 

chapter. This chopper or pulse width modulator is made up of two 


DC Power Supply 
Filter Capacitor 



Figure 3,1-1. Machine and Power Conditioner Schematic 






84 


transistor switches QM and Q6 as well as the two antiparallel diodes 
DB and DM respectively. 

The chopping action of these transistors (QM and QB) maintains an 
average d.c. current level in the attached coupling or chopper inductor. 
The time averaged valued of this current approximately equals the 
commanded machine current, IMG. This current command is directly propor- 
tional to the machine torque command because of the linear relationship 
between machine current and torque. The actual chopper inductor current 
does not track the command, IMG, exactly since this would require a 
chopper switching frequency of infinity. Instead, the chopper current 
oscillates in a sawtooth fashion of magnitude + ITOL amperes about the 
commanded current IMG. The antiparallel diodes, DB and DM provide 
paths for the decaying inductor current which results from the cut off 
of transistors QM or QB respectively. 

The modulated inductor current is then injected into the proper 
machine phase windings by a six leg transistor inverter bridge. Each 
of the three machine phase windings is serviced by two oppositely 
oriented transistor switches which allow control of the current 
entering or leaving a particular phase. In addition, each inverter 
transistor is spanned by a diode with a current orientation opposite 
to that of the attached transistor. These antiparallel diodes provide 
a path for ^he decaying phase currents which result when a phase is 
commutated (or disconnected) from the source. This greatly reduces the 
magnitudes of the voltage transients across the inverter transistors 
because of the lower time rate of change of the inductor current. 

The phase assignment of these switches is given below: 


85 


Phase A: Ql-Dl and Q4-D4 

Phase B: Q2-D2 and Q5-D5 

Phase C: Q3-D3 and Q6-D6 

Notice that the upper Inverter transistors Ql, Q2» and Q3 control 

positive phase currents (l.e. currents entering the machine), while 

the lower Inverter switches Q4, Q5, and Q6 control negative phase currents. 

The above mentioned commutation process results In a current unbalance 
at the node which connects the chopper Inductor and the top of the 
Inverter bridge. The commutation recovery diode, DR, provides a path 
for this excess current and at the same time clamps the maximum potential 
at this node to a level slightly above the supply potential. There 
would be a large voltage transient across the Inverter bridge If such 
a path did not exist. This Is a direct consequence of the large (L dl/dt) 
which results from the high resistance path. 

The logic which controls the switching of the chopper transistors 
QM and QB, and the logic which controls the six Inverter transistors 
Q1-Q6 Is described In detail In the chapter which follows. Mathematical 
representations of these components for the power conditioner and 
machine network model are developed In the next section. 


86 


3.2 MODELS OF THF, VARIOUS POWER CONDITIONER COMPONENTS 

This section deals with the derivation of equivalent electrical 
models of the various components of the power conditioner network. 

These components can all be classified as one of the five basic component 
types listed below: 

a) Battery 

b) Inductor 

c) Capacitor 

d) Switching Transistor 

e) Switching Diode 

The development of mathematical representations for these components is 
the objective of this section. 

The first component listed above is a battery which in the present 
case consists of a bank of 24 standard lead acid batteries connected 

■a 

together to provide a potential of approximately 270 volts. The actual 
discharge characteristics of such batteries are extremely difficult to 
model analytically. However at the present time such a detailed battery 
model is not warranted for two reasons. First the battery bank is 
usually operated at energy levels which are well below the levels where 
these discharge nonlinearities are apparent. And secondly this battery 
bank is of a temporary nature and therefore bears little resemblance 
to the actual d.c. supply present in the space shuttle. Based upon these 
assumptions the supply battery can be represented by a constant voltage 
source, Vg^^, in series with a resistance The constant voltage 

source, is assigned the value of the open circuit or no load 

battery voltage. The series resistance accounts for the internal 
battery resistance plus the resistance associated with the interconnec” 
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tlon between the battery and the power conditioner. A symbolic repre- 
sentation of this battery model Is given below: 

V$ 

• — f — • 

IDEAL BATTERY 



The relationship between Vg, and is given by ohms law 

as follows: 


' SatSaT '^BAT (3.2-1) 

Notice that the relationship between the assigned branch voltage and 
branch current obeys the consumer system of notation. Consequently, 
the battery receives or consumes power when the product ''^bat'^BAT 
positive, and likewise it supplies power to the load when ''^bat'^BAT 
negative. Henceforth, all branch current and voltage orientations 
will follow this convention. 

The power conditioner and machine model developed in this chapter 


contains a total of four inductances. The power conditioner network 
contains one of these four inductances, namely the inductance associated 
with the chopper inductor. This inductance is heavily influenced by 
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the magnitude of the chopper current. Figure (3.2-1) displays this 
nonlinear relationship between the chopper current and the chopper In- 
ductance. The computer program developed under this contract utilizes a 
cubic spline curve fitting algorithm to represent this nonlinear lity 
numerically. The remaining three Inductances are from the electrical 
model of the machine derived earlier In Chapter 2. These three leakage 
Inductances account for Che flux leakage of the stator phase windings. 

The volt-ampere relationship for the inductors is given by 

V - (3.2-2) 


where 

V is the voltage across the Inductor in volts, 

1 is the current through the inductor in amperes, 
L is the Inductance in henries and 
t is the time in seconds. 


A schematic representation of such an inductor branch is given below: 

L. 



A proper representation of the chopper and leakage inductors must 
include the voltage drop due to the non zero resistances associated 
with their windings. The addition of this resistive voltage drop results 
in the following; 




.V 



R is the winding resistance in ohms, 
w 


where 
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L is the winding resistance In henries, 
w 

Is the voltage drop In volts due to the winding 

Inductance L , 
w 

v^ Is the voltpge drop In volts due to the winding 

resistance R , 
w’ 

v Is the total voltage drop across the Inductor terminals 
In volts. 


This schematic can be expressed analytically as 


V ■ v_ + v„ 
L R 


V - L It- + R 1 (3.2-3: 

w dt w 

The power conditioner network as represented In this report contains 
only one capacitor. The functional relationship between the voltage and 
current of this capacitor Is given by: 

1 - C ^ (3.2-4: 

where 1 Is the current In amperes through the capacitor, 

C Is the capacitance In Farads, 

and V Is the voltage in volts across the capacitor terminals. 
Schematically this equation is represented by 


The transistors and diodes which make up the power conditioner 


basically function as simple ON/OFF switches. Consequently, their 




behavior can be approximated with relative accuracy by logical variables 
which take on values of either 1 (TRUE or ON) or 0 (’"'ALSE oi OFF) . The 
transistors are able to switch in very little time; typically in 
less than 1 us. For this reason the transistor and diode switches are 
represented by an ideal switch in series with a resistor to represent 
the voltage drop. The steps in the determiuation of the switch status and 
resistance are outlined below: 

STEP 1 . Determine the status of all power conditioner transistors. 
The transistor switching commands are generated by the pulse width 
modulator (chopper controller) and by the inverter controller (see chapter 
4.0). The commands are defined as follows: 

QMON Logical Switching Command for QM 
QBON Logical Switching Command for QB 
Q(l) Logical Switching Command for Q1 
Q(2) Logical Switching Command for Q2 
Q(3) Logical Switching Command for Q3 
Q(4) Logical Switching Command for Q4 
Q(5) Logical Switching Command for Q5 
Q(6) Logical Switching Command for Q6 

Each of these logical variables is either 1 (TRUE) or 0 (FALSE) corres- 
ponding to an ON or OFF switching command. 

STEP 2 . From the transistor switching commands generated in STEP 1, 
assign the proper resistance to each switch. Transistors commanded to 
be ON are assigned a low resistance based upon the voltage drop at 
rated currant. Transistors commanded to be OFF are assigned resistances 
which are 10,000 times larger than the ON resistance. 
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The procedure for the determination of the status of the switching 
diodes is analogous to the one used for the transistors. This procedure 
can be summarized as follows: 

STEP 1 . Determine the diode switching commands D(l), . . . , D(9) 
from the signs of the diode branch voltages. These commands are defined 
as follows: 

D(l) Logical Switching Command for D1 
D(2) Logical Switching Command for D2 
D(3) Logical Slwtching Command for D3 
D(4) Logical Switching Command for D4 
D(5) Logical Switching Command for D5 
D(6) Logical Switching Command for D6 
D(7) Logical Switching Command for D8 
D(8) Logical Switching Command for DM 
D(9) Logical Switching Command for DR 

Each of these logical variables is either 1 (TRUE) or 0 (FALSE) corresponding 
to the ON or OFF diode states. 

STEP 2 . Assign the proper branch resistance to each diode based 
upon the diode switching commands D(l), . . . , D(9) described in STEP 1. 
Diodes commanded to be ON are assigned a low resistance based upon the 
voltage drop at rated current. Diodes commanded to be OFF are assigned 
resistances 10,000 times larger than the corresponding ON resistance. 

The diodes which are connected across transistors (D1 D6, 

DB, and DM) are always OFF when the transistor is on because of their 
antiparalle connection. Consequently the transistor switching commands 
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take precedence over the corresponding antiparallel diode switching 
commands . 

The incorporation of these component itodels into a unified state 
space representation of the power conditioner and machine network is 
accomplished in the next section. 
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3.3 DERIVATION OF THE STATE EQUATIONS FOR THE POWER COITOITIONER 
AND MACHINE NETWORK 


3.3.1 Introduction and Modeling Approach 

The state equations of the combined power conditioner and machine 
network are derived in this section. These equations will be derived 
in standard state space form as symbolized by: 


• A *x + B *u 
-n -n —n -t 


where 


(3.3-1) 


and 


X is the vector of the network state variables, 

-n 

u is the vector of the network input or forcing functions, 

X is the vector cf the first derivatives of the elements 
— n 

of X with respect to time, 

A is the nonlinear coefficient matrix of the vector x , 

— n ~n 

B is the coefficient matrix of the vector u . 


The procedure selected for the derivation of these state equations 
is described in detail in References (1) and (5) . This procedure can 
be siusmarized briefly as follows : 

STEP 1 . Select a tree and cotree to represent the network topology. 
The rules for the selection of these is given in References (1) and (5). 
STEP 2. Derive the state equations in the intermediate form 


F X * G X + H u (3.3-2) 

"‘Tl "Tl Cl II 


F is the non singular coefficient matrix of x , 
-ici ^ 



where 
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G Is the nonlinear coefficient matrix of x , 
-n -n 

and 


H Is the coefficient matrix of u . 

-n — n 

STEP 3 . Convert the intermediate form of the state equations given 

in STEF 2 into the standard state space form of equation (3.3-1). This 

is accomplished by premultiplying both sides of equation (3.3-2) by 

F which results in 


X - [F • G ) X + [F • H ) u 
-n -n — n —n — n — n -n 

The matrices A and B of equation (3.3-1) can now be defined as 


-n -n 


-1 


G 

—n 


-1 


H 


(3.3-3) 


(3.3-4) 

(3.3-5) 


All derivations performed in these steps are kept in symbolic form, 
so as to reduce the time required for the numerical calculation of these 
equations . 
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3.3.2 Derivation of the Nonllaear Network State Equations 

Substitution of the component models developed in Section 3.2 into 
the power conditioner and machine schematic given In Figure (3.1-1) 
results In the network shown In Figure (3.3-1). The heavy lined branches 
In this figure represent the twigs of the chosen tree. The thin lined 
branches belong to the cotree, see References (1) and (5). The assumed 
direction of the curr'^'t In each of these branches is indicated by a 
solid arrow for the twxgs, and by the outline of an arrow for the links. 
The sense of the branch voltages with respect to the assumed current 
orientations Is determined by employing the consumer (load) system of 
notation. 

For purposes of Identification, each branch Is characterized by 
two different labels; a branch label and a branch component label. The 
branch label Identifies the branch number and Is denoted by the prefix B 
and the branch number. For example, B1 denotes branch one. The component 
labels are Identified by the prefixes C, E, L, and R (capacitor, emf, 
inductor, and resistor) and the component number. In addition to the 
branch labels, each node is identified by the prefix N and the chosen 
node number. Node NO is taken to be at ground or zero volts potential 
and is therefore used as the reference for the other node voltages. 

Branch currents are denoted by the prefix IB followed by the appro- 
priate branch number. For example, IB16 represents the current through 
branch 16. Similarly, the branch voltage is symbolized by VB16. 

The state variables for this network are the branch voltages of the 
twig capacitors and the currents through the link inductors. This 
corresponds tc VB5, IB21, IB22, and IB23, see Reference (1). 


g.RaXyi^ BIKR6 



Figure 3.3 1. The Machine and Power Conditioner Network Model 
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The steps in the derivation of the power conditioner and oachine 
state equations are listed below. 

STEP 1 . Determine the fundamental current loops for the given tree. 
The fundamental current loops are formed by taking one link at a time 
and tracing out the unique loop which would be formed by that particular 
link and the tree with all other links removed from the network. There 
are eight links in this network, consequently the number of fundamental 
current loops equals eight, see Figures (3.3-2), . . . , (3.3-9). 

STEP 2 . Determine the resistive link currents, IB16, . . . , IB20, 
in terms of the state variables, VB5, IB21, IB22, and IB23, and the 
sources, El, E2, E3, and E4. The link current, IB16, corresponding to 
link 1 (branch 16) is obtained by applying Kirchoffs voltage law around the 
loop through the tree formed by link 1. Solving for IB16 yields: 

IB16 » (VB5-E4)/R10 (3.3-6) 

Simllarily one obtains for link 2: 

IB17 = (VB5 + IB21R1)/(R1 + Rll) (3.3-7) 

The loop equations for links 3, 4, and 5 yield a set of coupled simultaneous 
equations which relate the link currents IB18, IB19, and IB20 to the 
state variables as follows: 

K • col {IB18, 1B19, IB20} *! L* 

^ -n 

col {-VB5, 1B21, IB22, 1B23} 
or in expanded form: 


R2 + R4 + R12 

R2 

R2 


IB18 

R2 

R2 + R5 + R13 

R2 

• 

IB19 

R2 

R2 

R2 + R6 + R14 ' 


IB20 1 












B7. y\2 



. Fundamental Current Loop Due to Link 1 (B16) 





















X A » PH,RS B3,E3 



Figure 3.3-3. Fundamental Current Loop Due to Link 2 (B17) 




Figure 3.3-4. Fundamental Current Loop Due to Link 3 (B18) 




Figure 3.3-5. Fundamental Current Loop Due to Link 4 (B19) 




Figure 3.3-6. Fundamental Current Loop Due to Link 5 (B20) 




Figure 3.3-7. Fundamental Current Loop Due to Link 6 (B21) 
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1 

R2 

R4 

R4 

1 

R2 

-R5 

0 

1 

R2 

0 

-R6 


VB5 


IB21 


IB22 


IB23 


(3.3-8) 


Premultiplying both sides of equation (3.3-6) 



yields : 


col {IB18, IB19, IB20) - C col {VB5, IB21, IB22, IB23} 

“Tl 

where C ■ k”^ • L 

"ti -n 

or In matrix form 


IB18 

IB19 

IB20 


Cll 

C12 

C13 

C14 

C21 

C22 

C23 

C24 

C31 

C32 

C33 

C34 


VB5 

IB21 

1B22 

IB23 


(3.3-9) 


STEP 3 . Form the loop equations corresponding to the link inductors. 
These loop equacions will involve the resistive link currents, IB16, . . ., 
IB21. These currents have already been expressed as a linear combination 
of the state variables ind the independent voltage sources; consequently, 
the loop equations for the link inductors can be written entirely in 
terms of these state variables and voltage sources. 

The loop equations corresponding to links 6, 7, and 8 (branches 21, 

22, and 23) become; 


Link 6 


L2 


dIB21 

dt 


VB5 • [ 


RI 

R1 + Rll 


+ R2(C11 + C21 + C31) - 1] 
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Link 7 


+ IB21 t-(Rl + R2 + R3) + — ——j + R2(C12 + C22 + C32) ] 

+ IB22 [R2(C13 + T23 + C33)] 

+ IB23[R2(C14 + C24 + C34)] (3.3-10) 


ai + L3)i2|2 + u^- 


dt 


VB5 ten R4 - C21R5] + IB21 [C21 R4 - C22RS] 


+ IB22 [C13 R4 - C23 R5 - (R4 + R5 + R7 + R8) ] 


+ IB23 t-(R4 + R7) + C14 R4 - C24 R5] 


+ El - E2 


(3.3-11) 


Link 8 

,, dIB22 . dIB23 _ 

Li + (LI + L4) - 

VB5 ten R4 - C31 R6] + IB21 [C12 R4 - C32 R6] 

+ IB22 [C13 R4 - C33 R6 - (R4 + R7)] 

+ IB23 [C14 R4 - C34 R6 - (R4 + R6 + R7 + R9)] 

+ El - E3 (3.3-12) 

STEP 4 . Write the capacitor current in terms of the state variables 
and independent voltage sources. Applying Kirchoff's Current Lav (KCL) 
at node 2 yields: 

Cl + IB18 + IB19 + IB20 
at 


- IB21 + IB17 + IB16 » 0 


(3.3-13) 
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Substituting for IB16, . 


Rewriting equations 


. . , IB21 in equation (3.3-12) results in: 

•'^5 'sfe ♦ STTm * 

•IB21 + C12 + C22 + C32 - 1] 

•IB22 [C13 + C23 + C33] 

•IB23 [C14 + C24 + C34] 

• E4/R10 (3.3-14) 

(3.3-10), (3.3-11), (3.3-12), and (3.3-14) in 


matrix form results in: 
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STEP 5 . The matrix equation given in equation (3.3-15) can be 
written symbolically as in equation (3.3-2) which is 


F X 

-n -n 


G X + H u 
-n -n 


This intermediate form of the network state equations can be transformed 
into the standard form given in equation (3.3-1) by premultiplying both 
sides of equation (3.3-15) by F where F ^ is defined below. 


















3.4 DERIVATION OF THE NETWORK OUTPUT EQUATIONS 


The network output equetlons relate the states and Inputs of the 
system to the branch voltages and currents. Knowledge of the four state 
variables and Inputs Is sufficient to uniquely determine the twenty-three 
branch currents and voltages. The algorithm for the calculation of 
these branch voltages and currents from the vectors and u^ Is outlined 
In the steps listed below. 

STEP 1 . Determine the first derivatives with respect to time of the 
state variables. These derivatives, k , are evaluated by means of 
equation (3.3-1). The vector x in matrix form can be written as: 



VB5 

IB21 

IB22 

IB23 


(3.4-1) 


STEP 2 . Determine the current IBS which flows through the twig 
capacitor Cl and the voltages VB21, VB22, and VB23 across the link 
inductors L2, L3, and L4 respectively, using ^ defined above as follows: 


IBS 


VB21 


VB22 


VB23 


Cl 

0 

0 

0 

0 

L2 

0 

0 

0 

0 

L3 

0 

0 

0 

_2J 

L4 


_d^ 

dt 


VBS 


IB21 


IB22 


IB23 


(3.4-2) 


STEP 3 . Calculate the voltage, VBIS, across the twig inductor LI as 
follows : 

VBIS * -LI • ~ (IB22 + IB23) (3.4-3) 

at 

STEP 4 . Calculate the currents through the resistive links (IB16, 

. . . , IB20) . From equations (3.3-6), (3.3-7), and (3.3-8), one obtains 
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STEP 5 . Calculate the twig currents IBl ZB15 from the 

link currents IB16, ...» ZB23 which were determined in the previous 
steps. These calculations make use of the network cut set matrix 
£ which is described in detail in References (1) and (5). The cue set 
matrix relates the twig current vector, 1^, and the link current vector, 
1., at the cut sets formed from the twigs. The twig current vector ^ 

“‘it “X 

consists of the twig currents IBl, . . . , IB15. The link current 
vector is composed of the link currents IB16, . . . , ZB23. The cut 
set matrix, £, can be partitioned as follows: 

where equals the identity matrix, of dimension 15, 

and is a 15 X 8 matrix as defined in Figure (3.4-1), Reference (5). 

Since ■ I^, equation (3.4-5) becomes: 

1, + • it - 0 (3.4-6) 

The tv±g currents ^ can now be written as 

it - • it (3.W) 

STEP 6 . Calculate the branch voltages across all branch resistances 

using the twig and link currents determined in the previous steps . The 

voltages VB6, . . . , VB14 across the twig resistors; Rl, . . . , R9, 
are given by 


iL 


- 0 


(3.4-5) 









,9 




VBIO 

VBll 

VB12 

VB13 

VB14 



IBIO 

IBll 

IB12 

IB13 

IB14 


(3.4-8) 


The voltages VB16, .... VB20 across the link resistors RIO, . . . , 
R14 are related to the link currents as follows: 


VB16 

VB17 

VB18 

VB19 

VB20 


RIO 

0 

0 

0 

0 

0 

Rll 

0 

0 

0 

0 

0 

R12 

0 

0 

0 

0 

0 

R13 

0 

0 

0 

0 

0 

R14 


IB16 

IB17 

IB18 

IB19 

IB20 


(3.4-9) 


Once all of these branch currents and '^oltages have been determined, 
one can calculate the Instantaneous power flowing through each of these 
branches by obtaining the product of the corresponding branch currents 
and voltages. The branch voltages are alsc used to update the status 
of the switching diodes. The next section deals with the instantaneous 
solution of the equations derived in this and the previous sections. 
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3.5 DISCRETIZATION AND SOLUTION OF THE CONTINUOUS TIME NONLINEAR 
NETWORK 

The nonlinear system of differential equations for the power 
conditioner and machine model must be Integrated forward in time to 
generate the time sequence of state vectors which characterize the 
dynamics of the EMA in time. The network equations are of the general 
form: 


Ax + Bu 


(3.5-1) 


This class of matrix equations can be numerically Integrated 
forward in time by considering time Intervals or steps, T, which are 
sufficiently small so as to allow holding A and ^ constant over such 
intervals, Figure (3.5-1). Given x(k) and u(k) at a time t - t^ - kx. 
one can calculate the state translstlon matrices £(‘T|^) and @('^j^) which 
relate s(t^^j^j to x(tj^) and u(k) as follows: 

x(tjj,^l) « *(xj^)x(tj^) + 0(Xj^)u(tj^) (3,5-2) 

The transition matrices 5 (Tj^) and §(tj^) are defined in terms of A(Tj^) 


and B(xj^) as follows: 


$(T ) s e-^ k « I + X, A(x) 


k 2 


(3.5-3) 


and 


rt»t 


S(V 


B dt 


t»t, 


or 


§< V ■ + 2^ + 3T * ooo)B(t^) 


where I is the identity matrix. 


(3.5-4) 
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The time steps Tq, *** are all chosen to equal the 
constant* t» In the EHA model. The steps followed by the computer 
program to Integrate the network forward by the N time steps are 
shown in the flow chart. Figure (3.5*2). This algorithm proceeds as 
follows: 

Step 1 . Determine the Initial system and transition matrices; 
-n^“^o^’ -n^V’ -n^V* -n^'^o^ respectively. 

Step 2 . Set the network integration counter, k, to zero. 

Step 3 . Increment the integration counter, k, by one. 

Step 4. Integrate the state variables forward in time from 
x^(tj^) to means of Equation (3.5-2). 

Step 5. Determine the network branch voltages, currents, and 
powers at time from the vectcTrs S^^^k+l^ -n^\+l^* 

Step 6. Determine if the network configuration has changed. 

NO: GO TO STEP 7 

YES: Update A„, B , « , and 0„. 

— n — n "u “n 

Step 7. Have there been N integrations? 

NO: GO TO STEP 3 

YES: STOP 
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4.0 STATE SPACE MODELS OF THE EMA SERVO 
LOOPS AND OF THE ROTATING MASSES 

The power conditioner and machine networks described and 
modeled in the previous chapter are controlled by four independent 
servo loops. These loops generate the control signals to switch the 
appropriate power conditioner transistors on or off depending upon 
the position command, the state of the machine line current, and the 
states of the rotating masses. 

The dynamics of this servo loop as well as the dynamics of the 
rotating masses are derived in this chapter. As before, these dyna- 
mics will be written in the standard state space form described earlier. 

The EMA to be modeled here is usually run with two channels 
braked. Reference (2). For this reason, the model developed is 
limited to one or two channel operation. The servo loop equations 
are highly nonlinear due to the current magnitude and race limiters 
which are used to control the magnitude and rate of the commanded 
machine torque. 

The first section of this chapter describes the configurations 
of the servo loops under one and two channel operation as well as 
Che configuration of the rotating masses. 

4.1 DESCRIPTIONS OF THE CONTROL LOOPS AND OF THE ROTATING MASSES 

A complete mechanization diagram of a single channel of the 
EMA is shown in Figure (4.1-1). Inspection of this diagram reveals 
the three basic feedback loops per channel. There is an outer posi- 
tion loop which feeds back the actuator position for comparison with 
the deflection command, DC. The time constant of the position error 


Velocity 



} 


Figure A.l-l. Mechanization of a Single EMA Channel 
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amplifier (AMPL D2<>7) is approximately 10 milliseconds. Inside this 
loop lies an inner velocity feedback loop which contains two amplifiers, 
labeled AMPL D8-1 and AMPL D2-10. Notice that the time constants of 
this loop are approximately one order of magnitude faster than those 
of the position error loop. The innermost loop consists of che current 
command magnitude and rate limiters, the power conditioner controller, 
the power conditioner, and the permanent magnet machine. This loop 
contains all of the nonlinearities accounted for in this effort. The 
time constants for this loop are in the microsecond range because of 
the rapid power conditioner dynamics. 

The wide variations in the time constants of these loops can leaa 
to numerical instabilities if one tries to solve the entire system of 
coupled differential equations in terms of a single matrix equation 
such as equation (4.1-1). Instead, the EMA will be broken into the 
three sub systems; 


X “Ax + B u 
-n -n-n -n-n 


X “Ax + B u 
-V -v-v -v-v 


(4.1-2) 

(4.1-3) 


and 

(4.1-4) 

“P “P“P “P”P 

based upon the magnitudes of the time constants. Equation (4.1-2) 
represents the dynamics of the power conditioner and machine network 
described earlier in Chapter 3. This portion of the EMA is outlined 
by block #1 in Figure (4.1-2). Equation (4.1-3) represents the por- 
tion of the EMA servo loop enclosed by block #2. This portion includes 
the velocity feedback loop as well as the nonlinearities introduced 
by the current command magnitude and rate limiters. These equations 



nr 
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J 


Figure A. 1-2. Classification of the EMA Components into Three Subsystems 
Based upon the Time Constants 
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are derived later in Section A. 3. The last system of differential 
equations symbolized by Equation (A. 1-4) models the dynamics of the 
rotating mass and of the position feedback loop, see block In 
Figure (A. 1-2). These equations are derived later in Section A. A. 

The servo loop gain constants KE and KR shown in Figure (A. 1-1) 
must be determined from tests performed upon the actual hardware at 
hand. This is accomplished by using the velocity gain and torque 
gain tests performed by Delco on the actual hardware, References 
(2), (3), (A). The results of these tests are presented in Figures 
(A. 1-3) and (A.l-A) respectively. 

These gain tests will be used to reformulate or consolidate the 
numerous gains of the velocity and position loops shown in Figure 
(4.1-2) into forms more suitable for the modeling process described 
later in this chapter. The proposed equivalent velocity and position 
loop structures are shown in Figure (A. 1-5). 

The gains KE and KV can be determined from the velocity gain 
test. Since the velocity loop can have only one gain, let KE ■ 1 
(dimensionless) for convenience. The other gain, KV, can now be 
determined from the velocity gain test. The velocity gain test was 
performed with the EMA output blocked and the position loop opened. 

The measured gain in the linear region was 0.035 amperes per rotor 
rpm, see Figure (A. 1-3). The input-output relationship of the velocity 
loop in the linear (undamped) region is given by 

= [KV • KEJRVEL (A. 1-5) 


From the gain test we have: 
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ICMDL , amps v 
RVEL ^rad/s«c‘' 


0 . 035 (^ • ) 


rpm 


radZsec' 


- 0.33422538 ( 


amps 

rad/sec"^ 


(4.1-6) 


Therefore; 


KV • KE - 0.33422538/5 

.066845076 (amps /rad /sec) (4.1-7) 

In a similar fashion, the gain, KERR, of the position error 

portion of the position feedback loop becomes: 

KERR • KE • 85. Kamps/deg flap) (4.1-8) 


The position feedback loop described above is lumped together 
with the rotating elements of the EMA. These elements are depicted 
in Figure (4.1-6). 

The next section describes and analyzes the EMA control logic. 
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4.2 EMA CONTROL LOGIC 

The EMA control logic which generates the switching cotamands: {QMON, 
Q60N, Q(l), Q(2), Q(3), Q(4), Q(5), and Q(6)} for the power conditioner 
transistors is described in this section. The aforementioned control 
logic can be divided into four basic components: 

a) The EMA Mode Decoding Logic 

b) The Current Command Magnitude and Rate Limiters 

c) The Pulse Width Modulator or Chopper Controller 

d) The Inverter Controller 

Figure (4.2-1) displays the functional interactions between these com- 
ponents. Descriptions and mathematical representations of these 
components follow immediately. 
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Figure 4.2-1. Black Diagram of the EMA Control Logic 
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4.2.1 EMA Mode Decoding Logic 

As stated earlier the Delco EM Actuator can operate In any one of 
the following inodes: 

a) The Motoring iw 'e 

b) The Regenerative Braking Mode 

c) The Plugging Mode 

These inodes of operation effect the operating configurations of both the 
control and power electronics as will be shown later in this chapter. 

The location of these inodes in the torque-speed plane was given 
earlier in Figure (2.2-1). Each of the six inodes shown in this figure 
is identified by a logical variable which is enclosed by a pair of 
parenthesis. These variables are defined as follows: 

MTRGl: Motoring - Speed (+) 

RGN4: Regenerative Braking - Speed (+) 

PLUG4: Plugging - Speed (+) 

MTRG3: Motoring - Speed (-) 

RGN2: Regenerative Braking - Speed (-) 

PLUG2: Plugging - Speed (-) 

Each of these variables is either TRUE or FALSE depending upon the operat- 
ing point in the torque-speed plane of Figure (2.2-1). Such a variable 
is TRUE when this point lies within a specified region and FALSE elsewhere. 
For example, MTRG 1 is TRUE whenever the motor torque and speed are both 
positive. 

These logical variables are generated by a decoder network from the 
rotor velocity RVEL (radians/secoud) and the commanded machine current 
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(or torque), ICMDl (aoperes). The decoder network is composed of two 
basic parts: 1) a comparator network which generates the logical 
variables {SPDHI, SPDLO, SPDPOS, SPDNEG, IPOS, and TNEG) from the 
rotor vellclty and the commanded current, and 2) a network of logic 
gates which operate on the outputs of the comparators to generate 
the mode variables {MTRGl, RGN4, PLUG4, MTRG3, RGN2, and PLUG 2}. 

Figure (4.2-2) shows the functional relationship between these two 
components of the EMA mode decoder network. 

The comparisons performed upon the analog signals, RVEL and ICMDl, 
by the comparator network are listed in the left hand column of Table 
(4.2-1). The corresponding logical output signals are given in the 
right hand column. 

The comparator output signals are then processed by the logic 
network shown in Figure (4.2-2). The outputs of this network are 
designated by the logical variables {KTRGl, RGN4, PLUG4, MTRG3, RGN2, and 
PLUG2} defined earlier. The boolean expressions, see Reference (2), 


1 define thl<: logic network are: 
MTRGl - TPOS • SPDPOS 

(4.2-1) 

RGN4 - TPOS * SPDHI 

(4.2-2) 

PLUG4 • TPOS • SPDPOS • SPDHI 

(4.2-3) 

MTRG3 » TNEG • SPDNEG 

(4.2-4) 

RGN2 = TPOS • SPDLO 

(4.2-5) 

PLUG2 » TPOS • SPDNEG • SPDLO 

(4.2-6) 


The dot ( • ) operator which appears in these equations denotes the 
logical AND operator. 



Figure A ,2-2. Functional Block Diagram of the 
Mode Decoder 










Table (4.2-1) Functional Description of the Comparator for the 
EMA Mode Decoder 


COMPARATOR FUNCTION 


COMPARATOR OUTPUT 

RVEL > VHI 

YES: 

SPDHI ■ 1 or True 


NO: 

SPDHI ■ 0 or False 

RVEL < VLO 

YES: 

SPDLO - 1 


NO: 

SPDLO - 0 

RVEL > 0 

YES: 

SPDPOS » 1 


NO: 

SPDPOS - 0 

RVEL < 0 

YES: 

SPDNEG « 1 


NO: 

SPDNEG ■ 0 

ICMDl > 0 

YES: 

TPOS • 1 


NO.- 

TPOS ■ 0 

ICMDl < 0 

YES: 

TNEG - 1 


NO: 

TNEG - 0 
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The six EMA mode status variables derived here later become inputs 
to the command current magnitude limiter, the chopper controller and 
the Inverter controller which are covered in the remainder of this 


section. 
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4.2.2 Current Coinmand Magnitude and Rate Limiter 

The permanent magnet brushless d.c. machines were designed to 
deliver maximum shaft torques of up to 13.558 nt-m (120 In- lb) each. 

This level of output torque corresponds to a current command, ICMDl, 
of approximately 60 amperes. This, unfortunately. Is equal to the 
current rating of the power conditioner switches. Consequently, at 
rated torque the switches are operating close to the maximum rated 
conditions, with the ever present danger of voltage and current 
transients of dangerous magnitudes. For this reason, the maximum 
machine torque during the motoring and regenerative braking modes has 
been limited to 9.0387 nt-m (80 In- lb). This torque limit Is accom- 
plished by clamping the current command, ICMDL, at +40 amperes. During 
plugging, ICMDL Is limited to 25 amperes because of the circulating 
currents which exist In the machine windings (see Section 2.2.3). 

These relationships between the undamped current command, ICMD, and 
the magnitude limited current command ICMDL, are Illustrated In Figure 
(4.2-3). 

Furthermore, to prevent rapid fluctuations in the torque or 
current command, the magnitude limited current command, ICMDL, Is passed 
through the rate limiter shown in Figure (4.2-4). The output of the 
rate limiter is designated by the variable ICMD. 

The rate limiter operates by Integrating the difference between 
the Input, ICMDL, and the output, ICMDl, Figure (4.2-4). Therefore, 
at steady state, the integrator input is zero which results in a clamped 
output, ICMDl, equal to ICMDL. 
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The maximum race of change of Che race limlced command currenc, 

XCMDl, can be calculaced by lecclng Che InpuC, ICMDL, cake on Che maximum 
values during each of Che chree modes. During Che mocorlng and regenera- 
clve braking modes XCMDL has a maximum value of 40 amperes. Assxmlng 
chac XCMDl Is inlclally zero, Che oucpuc of che clamp, x^, (see Figure 
4.2-4)) equals: 

^2 - Clamp (^) - 12 

Therefore che maximum race of change of XCMDl during mocorlng and regenera- 
Cion becomes: 

^ (XCMDl)j^ - 379 • 12 • 5 » 22740 Amps/Sec 
or 

•jr (ICMDDmav ■ 22.74 Amps/ms (4.2-7) 

at rliUv 

During Che plugging mode XCMDL is limlced to 25 amperes, therefore; 

X 2 “ Clamp (^) - 12 

and 

4 ~ (XCMDl) = 22.74 Amps/ms. (4.2-8) 

at 

The response envelope of the current command rate limiter is shown in 
Figure (4.2-5). The solid curves represent the maximum rate of 
change of XCMDl as calculated above and assuming open loop operation. 

The dashed lines represent expected outputs under closed loop operation. 
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The magnitude limited command current, ICMDL, and the rate and 
magnitude limited command current, ICMDl, control the operation of 
the chopper controller which is presented next. 


4.2.3 Current Chopper 


The EMA control loops generate a torque cosmand, TCMD, in response 
CO the position error, motor speed and the mode of operation. This 
torque command can be translated into an equivalent current command, 

ICMDl, in a linear fashion as follows: 

ICMDl - TCMD/0.1129 (4.2-9) 

This specified machine current is developed and maintained by a pulse 
width modulator or chopper network which is shown in Figure (4.2-6). 

The components of this network are: 

a) Transistor QM 

b) Transistor QB 

c) Diode DB 

d) Diode DM 

e) Coupling or Chopper Inductor L2 

Transistors QM and QB are used to maintain the line current, IM, within 
a specified tolerance of +IT0L (amperes) about the command current IMC. 

IMC is defined as: 

IMC » I ICMDl I . (4.2-10) 

The two diodes, DB and DM provide paths for Che decaying chopper 
inductor current which results when QB or QM respectively are switched off. 

The chopper network can take on two different operating configurations 
depending upon the EMA mode of operation. During Che motoring and plug- 
ging modes, QM is switched on and off while QB is always off. This 
allows current from the battery to be pulsed into the coupling (chopper) 
inductor . 

The dynamics of the active chopper elements during the motoring 
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and plugging modes of operation are displayed graphically In Flgurf 
(4.2-7). During these modes, the magnitude of the d.c. line current, 

III, flowing Into the Inverter Is monitored and compared with the crm- 
manded current, IMC. If IM Is below the commanded current, UlC, by 
a specified tolerance of ITOL amperes, QM Is switched on. This caust^s 
a rise Ir IM as shown In Figure (4.2-7). As soon as IM reaches a 
values (IMC + ITOL) amperes, QM Is switched off, which gives rise to 
a large voltage across the inductor terminals. This Induced voltage Is 
caused by the rapid time rate of change of the Inductor current which 
is produced by switching QM off. The rapidly rising inductor voltage 
turns on diode DM; thereby providing a low resistance path for the 
decaying inductor current. When IM decays to a value below (IMC- ITOL) 
amperes, QM is switched on again, thereby repeating the cycle. 

The regenerative braking case Is different from the motoring and 
plugging cases described above because the d.c. line current, IM, Is 
negative. In this case, transistor QB, Instead of QM as before, becomes 
the active element. QB switches on and off to maintain IM within the 
specified tolerance about the commanded current, ICMDL, Figure (4.2-8). 
IM is not rate limited as is the case during motoring and plugging 
because ICMDL instead of ICMDl is used as the current command. When 
QB is switched off, the decaying Inductor current flows through diode 
DB and into the supply battery; consequently, the kinectic energy of 
the rotating masses is converted into electrical energy. 




Control of these transistors Is accomplished by analog comparators 
and logic gates which are located In the low level electronics. A 
black box representation of this chopper controller Is shown In Figure 
(4.2-9). The first set of controller Inputs: {IM, IMC and ICMDL} are 
analyzed by two analog to digital comparators. The comparator outputs 
are designated by the logical variables QMC and Q6C. The Input-output 
relationships for these comparators are given In Figure (4.2-10). 

Notice that both comparators have a built In hysteresis to prevent 
toggling on Input noise. The operating locus defined by such hysteresis 
loops Is shown In Figure (4.2-11). The functional relationships 
between the comparator Inputs and outputs are given by 


QMC » 


and 


QBC = 


1 IF IM < IMC 
0 IF IM > IMC 


1 IF IM > -IMC 
0 IF IM < -IMC. 


(4.2-10) 


(4.2-11) 


The second set of Inputs {RGN2, RGN4} to the chopper control logic 
are processed together with QMC and QBC to generate the chopper switching 
commands, QMON and QBON, as follows: 


RGN » RGN2 + RGF4 

(4.2-12) 

QMON « QMC • ^ 

(4.2-13) 

QBON » QBC • RGN 

(4.3-14) 


These commands control the switching of the chopper transistors QM and 
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QB. The switching commands for the remaining power conditioner transis- 
tors {Ql, . . . , Q6} of the inverter bridge are obtained next. 
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4.2.4 Inverter Control 

Control of the six Inverter transistors {Ql, . . . , Q6} Is accom- 
plished by six logical, switching commands, donated by Q(l) through 
Q(6). The switching commands are generated from the RPS logic signals 
(AA, BB, and CC) and from the logic signals FF and RGN. FF represents 
the sign of the commanded torque (current) . FF is TRUE for ICMDl positive 
and FALSE for ICMDl negative. The logical variable, RGN, defined in 
equation (4.2-1) indicates whether or not the machine is functioning as 
a regenerative brake. 

The boolean expressions which generate Q(l) through Q(6) from AA, 


BB, CC, FF, 

and RGN, References 

(1) and (2), are given below: 


QAPl ■ 

AA * 

BB • FF 

+ AA • 

CC • FF 

(4.2-12) 

QANl » 

AA • 

BB • FF 

+ AA • 

W • FF 

(4.2-13) 

QBPl - 

BB • 

^ • FF 

+ AA • 

BB * FF 

(4.2-14) 

QBNl - 

BB • 

CC • FF 

+ AA • 

BB • FF 

(4.2-15) 

QCPl - 

QAPl 

• QBn 



(4.2-16) 

QCNl - 

QANl 

• QBNl 



(4.2-17) 

and finally 






Q(l) • 

QAPl 

• 



(4.2-18) 

Q(2) - 

QBPl 

• RGN 



(4.2-19) 

Q(3) - 

QCPl 

• S 



(4.2-20) 

Q(4) - 

QANl 

• RGN 



(4.2-21) 

Q(5) - 

QBNl 

' RGN 



(4.2-22) 

Q(6) - 

QCNl 

• RGN 



(4.2-23) 

Notice that 

Q(l) 

through 

Q(6) are always OFF during regeneration. 

These 
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relationships between the Input and output signals of the Inverter 
controller during the motoring and plugging modes are clearly displayed 
in Figures (4.2-12) and (4.2-13) respectively. 



FF 

AA 

BB 

CC 

RGN 

QAPON 

QBPON 

QCPON 

QANON 

QBNON 

QCaJOli 


0 


F 
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4.3 Derivation of the Velocity Loop State Equations 

The velocity loop, which was Introduced earlier In Section 
4.1.4, Is modeled here by means of the same state space approach used 
to represent the dynamics of the machine and power conditioner net- 
work. The transfer function form of this loop Is shown In Figure 
(4.3-1). The state equations developed In this section are based 
upon this figure. 

The transfer function representation of the velocity loop can 
be transformed into the signal flow graph of Figure (4.3-2). The 
outputs of the Integrators (the 1/S branches) are denoted by the 
state variables: XV(1), XV(2), and XV(3). The remaining state 

variables, XV(4) and XV(5), are Induced state variables at the out- 
puts of the clamps, and C£ 2 * functions performed by these 

clamps are symbolized by f^(*) f2(0» respectively. A descrip- 

tion of these clamps has been given earlier, in Section 4.2.2. 

The switches, SW^ and SW 2 at the inputs of the clamps, close 
at the beginning of each Integration interval. The data passed by 
these switches is immediately clamped for one entire Integration 
time step by C£^ and CJI 2 respectively. 

The differential equations which govern this network can be 
obtained by writing the equations for XV(1), XV(2), and XV(3) at 
Inputs to the three integrators as follows: 

jfvd) » KE • (PE - XV(3))/TAU1 - XV(1)/TAU1 (4.3-1) 

XV(2) = KRl • XV(5) (4.3-2) 

XV(3) = KV * RVEL/TAU2 - XV(3)/TAU2 (4.3-3) 
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XV(2)*'‘ ■ 

XV(2)" 

(4.3-8) 

XV(3)"'‘ - 

XV(3)“ 

(4.3-9) 

XV(4)'‘‘ - 

fj_[XV(l)"] 

(4.3-10) 

XV(5)‘^ - 

f2[XV(4)"] 

(4.3-11) 


In matrix form these equations become 




1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

M-) 

0 

0 

0 

0 

O 

0 

0 

CM 

0 



or symbolically 


x'^-Vx" (A. 3-12) 

The Interpretation of the state variables In terms of the velocity 
loop transfer functions Is given In Figure (4.3-1). 

The state equations for the mechanical and position loops are 
derived In the section which follows Immediately. 


4.4 DERIVATION OF THE STATE EQUATIONS FOR THE MECHANICAL AND POSITION 
LOOPS 

The mechanical and position loops consists of the first order 
flap position loop, depicted In Figure (4.1-5), and the rotating 
masses of Figure (4.1-6). 

The state equations for the rotating masses, with one or two 
active machines, can be obtained easily by representing this system 
by an equivalent mechanical network based upon the network elements 
of Figure (4.4-1). These elements are; 
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a) Torque, T (NT-M) 

b) Moment of Inertia, J (KG-M^) 

c) Stiffness Coefficient, K (NT-M/RAD) 

d) Viscous Damping Coefficient, B (NT-M/ RAD /SEC) 

The nodes of these branches correspond to the angular positions of the 
shafts, gears, rotors, etc. which make up the rotational system. 

Using these elements, the rotational system of Figure (4.1-6) 
is represented as shown in Figure (4.4-2). The angles 6^, 6^, 6^,, 

6q and 6p are defined in this figure. The change in the 
shaft torques through the different levels of the gearing is repre- 
sented by the gear ratios between the various shaft torques. 

The differential equations which describe this system are 
derived by summing the torques at the seven nodes designated by 

®M’ ®A* ®C’ ®1’ ®0’ * ®F’ ®F Figure (4.4-2). Performing 

these summations results in: 


TM - JM • + T1 

(4.4-1) 

T2 = T3 

(4.4-2) 

T4 * T5 

(4.4-3) 

T6 « T7 

(4.4-4) 

T8 - (6q - N3 • 0p)KACT 

(4.4-5) 

T9 - -(N3 • 6p - 6 q)KACT 

(4.4-6) 

TIO « JF • 6p + BF • 6p + KF • 6p 

(4.4-7) 


Also 


T2 


T1 • N1 


(4.4-8) 
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T4 

- T3 • NX 

(4.4-9) 

T6 

- T5 • N2 

(4.4-10) 

T8 

- 2T7 

(4.4-11) 

TIO 

- T9 • N3 

(4.4-12) 


Now the angle 6^ can be expressed in terms of 6^ as follows: 


2 • NX • N1 • N2 ®M 


(4.A-13) 


By means of equations (4.4-8) through (4.4-13), one can rewrite equa- 
tions (4.4-1) through (4.4-7) in terms of two second order differen- 
tial equations in the variables 6^ and 6^,, as follows: 

TM - JM • + TACT/2 • NX * N1 • N2 (4.4-14) 

N3 • TACT - JF • 6p + BF . 6^. + KF • 6j. (4.4-15) 

where 

TACT « T9 


These equations can be converted to the S-domain and therefore to 
a transfer function type block diagram. Taking the Laplace trans- 
forms of both sides of equations (4.4-14) and (4.4-15) and solving 
for S6^ and S6p yields: 



TM - TACT/2 • NX • N1 • N2 
SJM 


(4.4-16) 


S6p - N3 • TACT - BF • (S6p) - (KF/S) • (S6p) (4.4-17) 


TACT can be expressed in terms of S6j^ and S6p by the Laplace trans- 
formed equations (4.4-5) and (4.4-13) as follows: 
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TACT 


S6 



2 • NX • N1 • N2 


- N3 • (S6j.) 


KACT 

S 


(4.4-18) 


Combining equations (4.4-16) through (4.4-18) with the block diagram 
of the position loop In Figure (4.1-5) resultl fft*t!le func- 

tion block diagram of the mechanical and position loops. Figure (4.4-3). 

This block diagram can be readily translated Into the signal 
flow graph of Figure (4.4-4). The state variables of this network, 
s^nnbollzbii by XM(1) through XM(5), are the outputs of the Integrators. 
Writing the equations at the Inputs; 301(1) through 301(5), of the 
Integrators yields (In matrix form) : 


0 

-1 

4-Nl-N2*JM 

0 

0 

0 

KACT 

2-N1-N2-NX 

0 

-N3-KACT 

0 
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0 

N3 

JF 

BF 
" JF 
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1 
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KP-KERR 

TAU3 

1 

“tAU3 


301(1) 
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301(3) 

m 

301(4) 


301(5) 



XM(1) 


301(2) 


XM(3) 


XM(4) 


XM(5) 


1 

JM 

0 

0 

0 

0 

0 

0 

0 

0 

KERR 

TAU3 


TM 


DC I 


or symbolically as 


X 

P 


A X 

-p -p 


+ B u 
-P -P 


(4.4-19) 
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The procedure used to Integrate the EMA model forward in time 
is described next. 

4.5 SOLUTION OF THE OVERALL EMA SYSTEM STATE EQUATIONS 

The three basic subsystems of the EMA model developed in pre- 
vious sections and chapters are described by the matrix state equations: 

Equation (3.3-15) Power Conditioner and Machine Model 

Equation (4.3-16) Velocity Loop Model 

Equation (4.4-19) Mechanical and Position ?op Model 

The integration of the resulting fourteenth order state model forward 
in time requires a solution technique capable of handling nonlinear 
differential equations as well as interfacing the different integration 
speeds used for the three subsystems. The solution procedure which 
performs these tasks is the topic of this section. 

The signal paths between the various components of the EMA model 
as well as the component transfer functions (S-domain) are given in 
Figure (4.5-1). Notice that the power conditioner and machine models 
are not given in this diagram because of the nonlinear nature of these 
equations. The procedure used to integrate these equations forward in 
time is given in Section 3.5. 

The integration of the three EMA subsystem models is accomplished 
by one of two time discretization schemes. In the first schema, the 
three sets of differential equations are integrated forv;ard in time 
at the same rate or time step. This scheme is illustrated in Figure 
(4.5-2a). Notice that each block shown represents one time step. 

The integration order is indicated by the numbers within the blocks. 
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This scheme is inefficient and could result in excessive roundoff since 


the time constants of the three blocks are widely separated. To re" 
duce this inefficiency a second time discretization technique was 
adopted, Figure (4.5-2b). In this scheme, the time steps are scaled 
according to the relative speeds of the dynamics associated with the 
subsystems. Typical values of these time steps are: 


Power Conditioner and Machine Model: 5-20 vs 


Velocity Loop Model: 50-100 vs 


Mechanical and Position Loop Model: 100—500 vs 
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(b) Variable Time Steps 


Figure (4.5-2) Integration Options for the EMA Model 
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The staggering of the time steps, which is evident from the 
order of Integration (see Figure (4.5>2b), Introduces a problem con- 
cerning the generation of the Inputs for the slower Integration loops. 
That Is, some of the Inputs for the faster loops are functions of the 
state variables of the slower loops. These inputs are generated by 
Integrating the slowest loop first, and so on, and using linear Inter- 
polation to generate the proper Input sequence for the fast loops. 

The Interpolation scheme Is Illustrated by Figure (4.5-3). The 
solid points, (*), represent actual Integration points In time, while 
the points marked by an (x) are Interpolated. 

The number of integrations required for a complete step re- 
sponse can be greatly reduced replacing the detailed power conditioner 
and machine model of equation (3.3-15) by the following: 

TM - KT * ICMDl . (4.5-1) 

where TM Is the machine torque (Nt-m) 

KT is the machine torque constant (Nt-m/amp) 
and ICMDl Is the filtered current command (amps) 

In addition, the rotor velocity is clamped at YMSAT » +9500 rpm to 
account for the velocity saturation which occurs during motoring 
because of increasing phase back emfs. This simplification of the 
EMA model reduces the order of the state equations from fourteenth 
to tenth and at the same time it allows the use of large time steps. 

To further increase the program flexibility, an option has 
been added to allow steady state analysis of the power conditioner 
and machine network. Steady state, in this case, does not imply 
frequency domain analysis but rather it implies a constant torque 
command and a constant machine speed. These three program options 
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